Selection of the Ground State for Nonlinear 
Schrodinger Equations 



We prove for a class of nonlinear Schrodinger systems (NLS) having two 
nonlinear bound states that the (generic) large time behavior is character- 
ized by decay of the excited state, asymptotic approach to the nonlinear 
ground state and dispersive radiation. Our analysis elucidates the mecha- 
nism through which initial conditions which are very near the excited state 
branch evolve into a (nonlinear) ground state, a phenomenon known as 
ground state selection. Key steps in the analysis are the introduction of a 
particular linearization and the derivation of a normal form which reflects 
the dynamics on all time scales and yields, in particular, nonlinear Master 
equations. Then, a novel multiple time scale dynamic stability theory is 
developed. Consequently, we give a detailed description of the asymptotic 
behavior of the two bound state NLS for all small initial data. The methods 
are general and can be extended to treat NLS with more than two bound 
states and more general nonlinearities including those of Hartree-Fock type. 
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1. Introduction and statement of main results 

In this paper we study the detailed dynamics of the nonlinear Schrodinger equation 
with a potential (NLS) 1 : 

id t <f) = H(i) + \\<f)\ 2 <f). (1.1) 

Here, H = —A + V(x) is a self-adjoint operator on L 2 (IR 3 ) and A is a coupling 
parameter, assumed real and of order one. We assume that V(x) is a smooth 
potential, which decays sufficiently rapidly as \x\ tends to infinity (short range). 
Finally, we assume that the operator H has no zero energy resonance [23, 30], a 
condition which holds for generic V. 

NLS is a Hamiltonian system with conserved Hamiltonian energy functional: 

H en [0] = j\Vm\ 2j rV{x)\<i>{x)\ 2 + ±\<t>{x)\Ux (1.2) 

and additional conserved integral 

N[0] = J \<j>(x)\ 2 dx (1.3) 

These conserved integrals are continuous in the if 1 (]R 3 ) topology. An extensive 
discussion of the well-posedness theory can be found in the [6, 21, 47]. In partic- 
ular, NLS is well-posed globally in time in the space iJ 1 (!R 3 ), for initial data, 0o, 
which is sufficiently small in H l . 

Throughout this paper we shall assume that O has sufficiently small H 1 norm. 
We shall use the notation: 

£o = Uo\\hi (1-4) 

For V(x) which decays sufficiently rapidly as \x\ tends to infinity (short range 
potentials) the spectrum of H [34] consists of discrete spectrum, a<i{H), consist- 
ing of a finite number of negative point eigenvalues, and continuous spectrum, 
a c (H) = [0, oo). The dynamics of solutions for A = (linear Schrodinger) is very 

1 In a similar way, one can treat nonlincarities Aif [|</>| 2 ] is a nonlinear term which can be local 
or nonlocal. Typical examples are: if[|0| 2 ] = \<j>\ p and if[|(/>| 2 ] = /C*|</>| 2 , for some convolution 
kernel, K. 
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well understood. Let ipj* and Ej* denote bound states and bound state energies 
of the linear Schrodinger operator H: 

Htpj* = Ej^tpj*, {^j*,^k*)L 2 — 8jk- (1-5) 

Arbitrary initial conditions in an appropriate Hilbert space, evolve as t tends 
to infinity into a time-quasiperiodic part consisting of a superposition of time 
periodic and spatially localized states with frequencies given by the eigenvalues 
and a dispersive or radiative part, which decays to zero as t tends to infinity in 
appropriate spaces, e.g. L p , p > 2, L 2 ((x)~ a dx). 

In order to be more precise, introduce P c *, the orthogonal projection onto the 
continuous spectral subspace of H: 

P C J = f - E 0/W>^*- (1.6) 

3 

The solution of the linear Schrodinger equation can be expressed as 

e~ iHt (f>o = E (^*,0o> ^*e-^** + e-* Ht P c ^ . (1.7) 
j 

The time decay of the continuous spectral part of the solution can be expressed, 
under suitable smoothness, decay and genericity assumptions on V(x), in terms 
of local decay estimates [23, 30]: 

|| (x)- a e- iHt P c *(f>o ||l W < C<*>" § II Wfa IUw, (1.8) 
o~ > o"o > 0, and L 1 — L°° decay estimates [24, 57] 

|| e -<Ht P w 0o IU~(flS) < II 0o lUi(flS)- (1-9) 

For A 7^ the bound states of the linear problem persist, and bifurcate from 
the linear states at zero amplitude into branches of nonlinear bound states [35]. 
Of interest to us is the detailed dynamics of the nonlinear problem on short, 
intermediate and long time scales and, in particular, the manner in which the 
nonlinear bound states participate in the dynamics. In [35] variational methods 
were used to establish the existence and orbital Lyapunov stability of bound states 
which are local minimizers of H en subject to fixed N; see also [55, 15]. This result 
says that initial data which is close, modulo a phase adjustment, in H l to the 
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ground state remains H l close to a phase adjusted ground state for all time. The 
H 1 norm is closely related to the conserved Hamiltonian energy of the system and 
is insensitive to dispersive phenomena. Therefore, the detailed dynamics is not 
addressed by this result. For example, could the large time dynamics consist of a 
nonlinear ground state plus a small nonlinear excited state part? The main result 
of this paper implies that this cannot occur. 

For the nonlinear problem, the simplest question to consider is the case where 
H has only one simple eigenvalue and the norm of the solution is small. The 
detailed dynamics was studied in [38, 39, 32, 53]. Small norm initial data are shown 
to evolve into an asymptotic nonlinear ground state and a radiative decaying part. 

In this paper we study the multibound state variant of this question. We 
consider the specific case where H has two simple eigenvalues, _E * and E u . The 
linear Schrodinger equation then has two time-periodic solutions ip^e" 1 ^* 1 and 
■0>i*e~ iEl **, with Hipj* = Ej^ipj*, ipj* G L 2 . Therefore [35], NLS has two branches 
of nonlinear bound states bifurcating from the zero state at the eigenvalues of H, 
^ ao e~ iEot and ^ ai e~ iElt , with E L 2 satisfying 

/M/„ • AM/,, -vl/„ /^vl/„ . (l.io) 

Here, otj denotes a coordinate along the j th nonlinear bound state branch and 

E 3 = E^ + 0(\a 3 \ 2 ). (1.11) 

In contrast to the linear behavior (1.7) our main result is the following: 

Theorem 1.1. Consider NLS with V (x) a short range potential supporting two 
bound states as described above. Furthermore, assume that the linear Schrodinger 
operator, H, has no zero energy resonance. 

(1) Assume the initial data, 0(0), is small in the norm defined by: 

[0(O)]x = HW(0)|U, (1-12) 

where k > 2 and a > are sufficiently large. Let (pit) solve the initial value 
problem for NLS. 

Assume the (genetically satisfied) nonlinear Fermi golden rule resonance condition 2 
r w , = A 2 tt (Ml, 5(H - > (1.13) 

2 The operator / S(H — w*)/ projects / onto the generalized eigenfunction of H with 
generalized eigenvalue w*. The expression in (7.4) is finite by local decay estimates (1.8); see 
e.g. [42]. 
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holds, where 

u* = 2E u -E Ot >0. (1.14) 

Then, as t — > oo 

0( t ) e-^W* a . (oo) + e lA V + , (1.15) 
in L 2 , where either j = or j = 1 . The phase ouj satisfies 

Uj(t) = u°°t + C*(\ogt) (1.16) 

Here, ^^(oo) JS a nonlinear bound state (section 3), with frequency Ej (oo) near 
Ej*. When j = 0, the solution is asymptotic to a nonlinear ground state, while in 
the case j = 1 the solution is asymptotic to a nonlinear excited state. 
(2) More specifically, we have the following expansion of the solution 4>(t): 

0(t) = e-^W* ao(t) +e-^W* ai(t) 

+ me-^^PcMt) + R loc (t) + Rm oc (t), 

where as t — > oo 

Uj(t) - Uj(t) -> 0, ^(t) - aj(oo) -> (1.17) 
and such that for each initial state, 0(0), 

|a (oo)| • |c*i(oo)| = 0, 

4>+(t) -> 0+ in Z 2 , 

ll^o C (-,t)|| 2 = 0(H) 

||^ OC (-,t)||oo = 0(t--*). 

Here, Ho(oo) is a small spatially localized perturbation of the operator 

a 3 (-A + V(x)) 

and P c = P c (TLq(oo)) , the projection onto its continuous spectral part. Finally, 
7Ti maps the vector (zi,z 2 ) to (21, 0). 



Remark 1.1. (A) Theorem 1 . 1 implies the absence of small norm time - quasiperi- 
odic solutions for this class of nonlinear Schrodinger equations [36]. Intuitively, 
one can explain why one expects only a pure state in the limit t — > 00 and how the 
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condition on uj* = 2Eu — Eo* arises. Our intuition is based on viewing the nonlin- 
earity as linear time-dependent potential; see also [36, 37]. An approximate super- 
position of a nonlinear ground state and excited state ~ ty ao e~ lEot + ^/ ai e~ lElt 
can be viewed as defining a self-consistent time-dependent potential: 



for \ai | << |oK j - As shown in [43], [27], [28] data with initial conditions given 
by the unperturbed excited state decay exponentially on a time scale of order 
t ~ O(\a ai\~ 2 ) provided the forcing frequency, E\ — E ~ E u — E * > —E u or 
lu* = 2E\* — Eq* > 0. 

(B) Theorem 1.1 implies asymptotic stability and selection of the ground state 
for generic small data. Part 1 of Theorem 1.1 implies a form of asymptotic com- 
pleteness. 

(C) Since we control the decay of solutions in W k '°° , our results imply global 
existence of small solutions in H s for all s sufficiently large. 

(D) The asymptotic state where |ai(oo)| ^ (and therefore |a (oo)| = 0) is non- 
generic. This can be seen by linearization about the excited state. The linearized 
operator, Hi, is a localized perturbation of an operator having embedded eigen- 
values in its continuous spectrum, under our hypothesis cj* = 2Eu — Eo* > 0. The 
connection between embedded eigenvalues in the continuous spectrum of an ap- 
propriate linear operator and the non-persistence of localized time periodic states 
and between embedded eigenvalues in the continuous spectrum of an appropriate 
linear operator was explored first in [36, 37]. It is well known that embedded 
eigenvalues in the continuous spectrum are unstable to generic perturbations; see, 
for example, [7, 14, 42]. In this case, the embedded eigenvalues are perturbed to 
complex eigenvalues, with corresponding eigenstates whose evolution is exponen- 
tially growing with time, under the condition (7.4). The perturbation to the linear 
operator with embedded eigenvalues is however both nongeneric (in that it comes 
from linearization of a Hamiltonian nonlinear term about a critical point of the 
energy) and breaks self-adjointness with respect to the standard L 2 inner product. 
A second order perturbation theory calculation shows that if cj* = 2Eu — E * > 
generically the embedded eigenvalue perturbs to an exponential instability [40]. 



W(x,t) 




A |^ ao + ^ ai e-^- E ^\ 2 

\\^ ao \ 2 + 2\ COs((£i - E )t + 7)*|a |*|ai| 



(1.18) 
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This suggests the existence of an unstable manifold of solutions for the nonlinear 
equation. The existence of such non-generic solutions of NLS with £1(00) 7^ for 
the full nonlinear How has recently been demonstrated [50]. 
(E) Theorem 1.1 is stated for the case of two nonlinear bound state branches. 
The technique of proof, however, can be used to consider the more general case. 
We expect results which are analogous to those of our main theorem, but more 
complicated due to the presence of: direct bound state - bound state interactions, 
bound state - continuum interactions and bound state - bound state interactions 
mediated by the continuum. Multimode Hamiltonian systems have been consid- 
ered in the context of linear time almost periodic perturbations have been studied 
in [28]. 

Relation to other work 

We now wish to further put our results in context. Research on nonlinear scat- 
tering in the presence of bound states has followed two related lines. 
A. Nonlinear dispersive waves in systems with defects, potentials,... 

Our analysis centers around nonlinear bound states which bifurcate from linear 
bound states of the operator, H, obtained by linearizing about the zero solution. 
These bound states exist at all sufficiently small amplitudes (measured in any H s 
norm, s > 0). The behavior of the "bifurcation diagram" for larger amplitudes 
depends in a detailed way on the details of the nonlinearity, the spatial dimension 
and the norm [35]. Such nonlinear bound states are also called nonlinear defect 
modes, nonlinear localized modes or nonlinear pinned modes. They are localized 
about or "pinned" to the support of the potential, V, and arise due to a local de- 
viation from translation invariance or a "defect" in the homogeneous background 
which acts as an attractive potential well. To get a more refined picture of the 
dynamics than in the H 1 theory, one must consider the linearized evolution about 
the family of nonlinear ground states. This linearized operator has continuous and 
discrete spectral parts inherited from the linear bound state spectral structure. 
In particular, the discrete spectrum contains an eigenvalue at zero corresponding 
to the ground state and a pair of eigenvalues (located symmetrically about zero) 
corresponding to the excited state. Thus, at linear order a solution infinitesimally 
close to the ground state formally appears to be quasiperiodic in time - a ground 
state plus a small excited state oscillation. However, at higher order in perturba- 
tion theory one finds nonlinear resonant coupling of the neutral oscillatory modes 



8 



to the continuum and as a result these slowly damp to zero; generically, for very 
large time energy splits between the ground state and dispersive parts of the so- 
lution. This mechanism for relaxation to the ground state was earlier considered 
for the nonlinear Klein-Gordon wave equation with a potential, where the decay 
of "breather-like" solutions was studied [44]. In this work, small norm solutions 
relax to the zero solution via resonant energy transfer out of the bound state to 
radiation modes and dispersive radiation of energy to infinity; the zero solution 
plays the role of the ground state. Results concerning special classes of initial 
data are considered in work of S. Cuccagna [11] and T.-P. Tsai and H.-T. Yau 
[48, 49]. Similar results to those of Theorem 1.1 are considered in [51]. 
B. Nonlinear dispersive translation invariant equations: A closely related 
line of research focuses on the translation invariant nonlinear Schrodinger equa- 
tion. Here the equation is (1.1) with V taken to be identically zero nonlinear 
coupling parameter A < 0. In this case, the equation has solitary wave solutions, 
obtainable by minimization of 7i cn [0] subject to Af[(f>] = Ao- For NLS in dimen- 
sion n with cubic nonlinearity replaced by the general power nonlinearity 
we have that if p < 1 + -, the foregoing variational problem has a unique (up 
to translation) radially symmetric ground state solution for any Ao > 0. In the 
case when V(x) ^ is a potential supporting bound states, the small A/o solu- 
tions agree with the bifurcating bound states discussed above [35]. As pointed out 
earlier, constrained energy minimizers are H l orbitally Lyapunov stable [12, 55]. 
An interesting feature of solitary waves in the translation invariant case is the 
presence spurious neutral oscillations. These are sometimes called internal modes 
[22] . To explain this, consider the linearization about a ground state solitary wave 
(p < 1 + -). Due to the underlying symmetric group of the equation (transla- 
tion invariance, phase invariance, Galilean invariance,..) this linearization has a 
(generalized) zero eigenvalue of multiplicity related to the dimension of equation's 
symmetry group. In the non-integrable cases (n — 1, p ^ 3 and n > 2) the lin- 
earization has additional neutral modes. These neutral modes approach zero as p 
approaches 1 + - ; the dimension of the zero subspace jumps by two at p — 1 + - , 
the critical case, corresponding to the larger group of symmetries and the existence 
of a pseudo-conformal invariant [54]. V.S. Buslaev and G. Perel'man [3] consid- 
ered the problem in one space dimension and showed that nonlinear resonance 
of these "internal modes" with the continuum is responsible for their damping 
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on long time scales and the asymptotic stability of solitons. See also the recent 
work of V.S. Buslaev and C. Sulem [4]. Their analysis was restricted to one space 
dimension only in their use of explicit eigenfunction expansion methods to obtain 
the required local energy decay estimates. S. Cuccagna [10, 11] extended their 
results to more general nonlinearities and general space dimensions. In his anal- 
ysis, the required dispersive estimates are obtained by adapting K. Yajima's [57] 
approach in which the wave operators, which conjugate the linearized operator on 
its continuous spectral part to the constant coefficient "free" dispersive evolution, 
are shown to be bounded on W k,p spaces. This method was also used in [44]. 

Another feature, common to problems of type A. and B. is the use of the 
method of normal forms. In the context of nonlinear scattering, normal form 
ideas were used to obtain the local behavior in a neighborhood of a soliton in [3] 
and for the decaying breather-like state in [44]. In contrast to the normal form 
for finite dimensional Hamiltonian systems, resonant interaction with the contin- 
uous spectrum gives rise to a more general normal form which captures internal 
damping, due to energy transfer out of certain discrete modes to the continuum 
modes; see the discussion in the introduction to [44]. In the present work, we 
derive a nonlinear master equation, coupled equations for the renormalized (up to 
near identity transformations on the complex discrete mode amplitdues) discrete 
mode square amplitudes ("mode powers"), which governs generic dynamics on 
large intermediate and very long time scales. Normal forms of this type, expected 
to be valid for very long times, were derived and studied in the local analysis 
about the steady and " wobbling " kink-like solutions of discrete nonlinear wave 
equations in [25]. A key feature of the normal form of the current work is our 
analysis of its behavior on different time scales and the analysis of its transitional 
behavior across time scales, for general initial data. 

Finally, we point out that there are many important areas of application which 
motivate the for study of the class of models we treated in this paper. We mention 
two. At the most fundmanental level the nonlinear Schrodinger equation arises 
as a mean field limit model governing the interaction of a very large number of 
weakly interacting bosons [20], [46], [29], [9]. At a macroscopic level, it has been 
shown that equations of this type arise as the equation governing the evolution 
of the envelope of the electric field of a light pulse propagating in a medium with 
defects. See, for example, [17, 18, 19]. 
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2. Structure of the proof 

We now sketch our analysis. Certain notation is defined in the appendix of section 
16. In analogy with the approach introduced in [38], [39] in the one bound state 
case, we represent the solution in terms of the dynamics of the bound state part, 
described through the evolution of the collective coordinates a (t) and cxi(t), and a 
remainder <f> 2 , whose dynamics is controlled by a dispersive equation. In particular 
we have 

0(t, x ) = e -iSlMs)ds-m) ( ^ + $aiw + Mt ^ x) ) (21) 

We substitute (2.1) into NLS and use the nonlinear equations (1.10) for fy a . 
to simplify. Anticipating the decay of the excited state, we center the dynamics 
about the ground state. We therefore obtain for <3> 2 = (02,02) T the equation: 

id t ®2 = H (t)$ 2 + G(t, x, $ 2 ; d t a(t), d t a, d t Q(t)) (2.2) 

where, 7io(t) denotes the matrix operator which is the linearization about the 
time-dependent nonlinear ground state ^ ao (t)- The idea is that in order for fait, x) 
to decay dispersively to zero we must choose ao(t) and ai(t) to evolve in such a 
way as to remove all secular resonance terms from Q. Thus we require, 

Pc(Ho(t))$ 2 (t) = 0, (2.3) 

where P c (7io) denotes the continuous spectral projection of 7Yo; see also condition 
(5.13). Since the discrete subspace of Ti^t) is four dimensional (consisting of 
a generalized null space of dimension two plus two oscillating neutral modes) 
(2.3) is equivalent to four orthogonality conditions implying to four differential 
equations for ao, ci\ and their complex conjugates. These equations are coupled 
to the dispersive partial differential equation for $2- At this stage we have that 
NLS is equivalent to a dynamical system consisting of a finite dimensional part 
(6.23-6.24), governing dj = (ctj, aj), j = 0, 1, coupled to an infinite dimensional 
dispersive part governing $ 2 : 

— * 

id t a = A(t)a + F a 
id t $ 2 = W (*)$2 + i^ (2.4) 
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We expect A(t) and TCo(t) to have limits as t — > ±00. Our strategy is to fix 
T > arbitrarily large, and to study the dynamics on the interval [0, T\. In this 
we follow the strategy of [3, 11]. We shall rewrite (2.4) as: 



id t a = A(T)a + (A(t) - A(T))a + F a 
id t $ 2 = H(T)$ 2 + (n (t)-'Ho(T))$2 + F <l , (2.5) 

and implement a perturbative analysis about the time-independent reference lin- 
ear, respectively, matrix and differential, operators A(T) and 7io(T). 

More specifically, we analyze the dynamics of (2.5) by using (1) the eigenvalues 
of A(T) to calculate the key resonant terms and (2) together with the dispersive 
estimates of e^^ 1 P C (T) [10]. 

Note also that P c (T)$ 2 (t) ^ $ 2 (0 because $ 2 (i) e Range P c (i) ^ Range P C (T). 
We therefore decompose $ 2 as: 

$ 2 = disc(t; T) + r], V = P C (T) V (2.6) 

where disc(t; T) lies in the discrete spectral subspace of H.q{T) and show that 
disc(t; T) can be controlled in terms of rj. 

The expected generic behavior of this system is that ot\ and $ 2 decay with a 
rate t~*. This slow rate actually leads to an equation for a\ with the character 
dtOL\ ~ (7P — <9i@)ai + integrable in t; see (6.12). Thus, is chosen to satisfy 
d t Q ~ ensuring that a\ has a limit. In this way a logarithmic correction to 
the standard phase arises; see (1.16). 

Next we explicitly factor out the rapid oscillations from a,\ and show that, 
after a near identity change of variables (ao,«i) 1— > (a , that the modified 
ground and excited state amplitudes satisfy the system: 

i<9 t a = (C1022 + *r^) | 4 Qf + F a [a ^ u -q,t] 

i3t& = (cii2i-2ir w )|a riA| 2 A + ^[d 0) A, 77, t]; (2.7) 

see Proposition 7.1. It follows that a Nonlinear Master Equation governs Pj = 
Icijl 2 , the power in the j th mode: 

!U ' ] - 2FP?P • /.',,(/) 



: -ATPfPo + R^t). (2.8) 



d^i ,.,02 



dt 
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Coupling to the dispersive part, $2, is through the source terms Rq and R\. The 
expression "master equation" is used since the role played by (2.8) is analogous 
to the role of master equations in the quantum theory of open systems [13]. 

A novel multiscale Lyapunov argument is implemented in section 8 charac- 
terizing the behavior of the system (2.8) coupled to that of the dispersive part 
on short, intermediate and long time scales. We consider the system (2.8) on 
three time intervals: I = [0,t ] (initial phase) I\ = [t ? ^1] (embryonic phase) and 
I2 = [ti, 00) (selection of the ground state). 

For t > to, the terms Ro(t) and Ri(t) are shown (Proposition 12.1) to have 
the form 

Ro(t) ~ + Po(£o,t)PoP? (2.9) 

Rl(t) ~ + SJfyyfPtf? + Pl (£ ,t)P P?, (2.10) 

where 

b = O((t )- 1 ) , b 1 = O{{t )-$) (2.11) 

Ro(t) and Ri(t) have parts which are local in time and nonlocal in time. The 
handling of nonlocal terms is explained in section 8. 
We set (Proposition 9.2) 

Qo = Po ~ |y, Qi = Pi + |y (2.12) 

where b and b\ are positive and satisfy (2.11) as well. 

Discussion of time scales 

Estimates (2.9-2.10) and the definitions (2.12) imply an effective finite-dimensional 
reduction to a system of equations for the "effective mode powers": Qo(t) and 
Qi(t), whose character on different time scales dictates the full infinite dimen- 
sional dynamics, in a manner analogous to role of a center manifold reduction of 
a dissipative system [5]. 

Initial phase - t E I = [0, t ] : Here, J is the maximal interval on which Qo(t) < 0. 
If t = 00, then Po(t) = 0({t)~ 2 ) and the ground state decays to zero. In this 
case, we show in section 15 that the excited state amplitude has a limit as well 
(with may or may not be zero). This case is nongeneric. 
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Embryonic phase - t e h — [to,ti\: If to < oo, then for t > t : 



> 2r'Q Q2 

- < -WQoQj + O(Jq~ QT)- ™>4 (2.13) 



Therefore, Qo is montonically increasing; the ground state grows. Furthermore, if 
Qo is small relative to Qi, then 

j^- is montonically increasing, (2-14) 

Qi 

in fact exponentially increasing; the ground state grows rapidly relative to the 
excited state. 

Selection of the ground state t e h — [ti, 00): There exists a time t = t±, t < 
ti < 00 at which the 0(y/QoQ™) term in (2.13) is dominated by the leading 
("dissipative") term. For t > t 1 we have 

dQo 



dt 
dQ 



> 2T'Q Ql 
- < -^T'QoQl (2.15) 



dt 

It follows that Qo(t) — * Qo(oo) > and Qi(t) — > as t — > 00; t/ie ground state is 
selected. 

3. Linear and nonlinear bound states 

In this section we introduce bound states of the linear (A = 0) and nonlinear 
(A 7^ 0) Schrodinger equation (1.1). 

Bound states of the unperturbed problem 

Let H = —A + V(x). We assume that V(x) is smooth and sufficiently rapidly 
decaying, so that H defines a self-adjoint operator in L 2 . Additionally, we assume 
that the spectrum of H consists of continuous spectrum extending from to 
positive infinity and two discrete negative eigenvalues, each of multiplicity one. 

a(H) = {E *,E U } U [0,oo) (3.1) 

Therefore, there exist eigenstates, ipj* G V(H), j — 0, 1 such that 

Hipj* = Ej^ipj*. (3.2) 
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We also introduce spectral projections onto the discrete eigenstates and con- 
tinuous spectral part of H, respectively: 

p j*f = (^j*,f)^j*, J =0,1 

Pc* = I Pi 0* Pi* 

Nonlinear bound states 

We seek solutions of (1.1) of the form 

= e~ iE ^ Ej (3.3) 

Substitution into (1.1) yields 

HV Ej + A|* B . | 2 # B . = Ejijj Ej (3.4) 

We introduce a bifurcation parameter, otj, for the j th nonlinear bound state branch 
and define 

aj = {pij,aj) (3.5) 



Proposition 3.1. [35] For each j = 0, 1 we have a one parameter family, ty a . = 
tyj, of bound states depending on the complex parameter a j = |aj|e* 7j and defined 
for \otj\ sufficiently small: 

Vj(x) = a j i/} j (x;\a j \ 2 ) = e vl ' j \a j \^ j (x;\a j \ 2 ) 
= otj ( ipj*(x) + X\a j \ 2 ipj 1 \x; \aj\ 2 ) ) 

= aj ( if>j*(x) + X\aj\ 2 ^f\x] 0) + 0(X 2 \a j \ 4 ^j . (3.6) 



Here, 



^\x;0) = -(H - EtoT\l - P Ej .)tf. (3-7) 
Ej = Ejipij) = E^ + la^E^da^ 



= E^ + xJ^Jx \ aj \ 2 + 0(\aj\ 4 ). (3.8) 
The mapping ctj i— > (Ej(dj), ^(-; <%)) is smooth. 
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The proof uses standard bifurcation theory [31], which is based on the implicit 
function theorem. The analysis extends to the case of nonlocal nonlinearities. 
A variational approach can also be used to construct nonlinear bound states. 
Variational approaches, though more global, do not directly yield the information 
we require concerning smooth variation with respect to parameters. 

Remark 3.1. tyj depends on otj and otj. We shall compute derivatives of the 
nonlinear bound states \l/ a with respect to otj and Wj and use the notation: 

V j =(d aj ,d^) = (d j ,cT j ). (3.9) 

In what follows we shall "modulate" these bound states. That is, we shall 
allow a to vary with time. For convenience, we shall use the notation: 



*j{t,x) = #cy(t)(aO, 

Ej(t) = EjQajW 



4. Linearization about the ground state 

Let ^ denote a nonlinear bound state (ground state or excited state) of (1.1); see 
section 3. Then, 

+ A|^| 2 ^ = EV. (4.1) 
We first derive the linear stability problem. Let 

= + p)e~ lEt , (4.2) 

where p denotes the perturbation about Substituting (4.2) into (1.1) and 
neglecting all terms which are nonlinear in p and p, we obtain the linearized 
perturbation equation 

id tP = ( H-E + 2\\m\ 2 ) p+ \^ 2 p (4.3) 

Since p appears explicitly in (4.3) it is natural to consider the system for 

p = (;)■■ (^) 

id t p = Hp, (4.5) 
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where 

H = 4 A* 2 H-B + 2XW)- (46) 

where 

<*=( 1 n " ) ■ (4-7) 



Later in this paper we shall refer specifically to the linearization about a 
"curve" of bound states (\&j(t), Ej(t)) and will denote by Hj(t) the operator (4.6) 
with E replaced by Ej(t) and \& replaced by ^(t). Our main focus will be on the 
operator family 

H (t) = H Eo (t)^ (t)- (4.8) 

The nonlinear bound state \1> is linearly spectrally stable if the spectrum of 
H, cr(H), is a subset of the real line 3 . \I/ is linearly dynamically stable if, in an 
appropriate space, all solutions of the initial value problem for (4.5) are bounded 
in time. That is, in some norm e~' tHot is a bounded operator. Linear dynamical 
stability of the ground state follows from [54]. For this result and the necessary 
stronger dispersive estimates on e~ lHot [10], we require information on the discrete 
spectrum of 7i and the corresponding spectral subspaces. 

Before stating these results we observe that the operator e~ tHot can be ex- 
pressed in terms of the operator treated explicitly in [54], [10]. To see this, express 
the ground state as ^ = [^ol^ 7 Q^ol > 0), where 7 = arga is a constant. Set 
p = e^q. Then, by (4.3) we have: 

id t q = (H -E + 2\\V \ 2 ) q+ X^q (4.9) 
Now let q — u + iv, where u and v are real. Then, 



J Ho ! , (4.10) 






1 andiJo= tfn-£o-3A|*ol 2 (4 - U) 



3 The spectrum of H. has the symmetries one expects for Hamiltonian systems. The mappings 
A 1 ► —A and An A send points in the spectrum to points in the spectrum. Note that if £ is 
an eigenvector of H with eigenvalue \i then is an eigenvector of H with eigenvalue —fi. 
Therefore, t-^ = a^. 
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Note that 

pit) = mpit) = 7r ie - mot p 

= ( ir ie jflot Kg + in 2 e J&ot 3g ) , (4.12) 

where tti(z 1 ,z 2 ) = (z 1 ,0) and ^2(^1,^2) — (0, ^2)- Therefore, we have: 

Proposition 4.1. Estimates on e~ Mot are equivalent to those for e JHot and are 
independent of 7. 

We now turn to a detailed discussion of the spectral properties of Ho- 

Proposition 4.2. Consider Ho, the linearization about the ground state. Let a 
he sufficiently small. 

(1 ) a (Ho) is a subset of the real line. 

(2) cr discrctc (Ho) = {-/i, 0,/i}, where < p < \E \. 

(3) Zero is a generalized eigenvalue of Ho- The generalized null space, N 9 (Ho), is 
given by 

iV 9 ( Wo) = span (4 §£ ),(§|§£)1 (4,3) 



(4) ±/x are simple eigenvalues. We denote their corresponding eigenfunctions by 
£ M and For \cx \ small we have the expansion: 

p = E, - E + O(\a \ 2 ) (4.14) 



= (4-16) 

where c\(a) and 02(a) are real analytic functions in a. 
(5) cr(H ) - cr discTCtc (H ) = (-00, E ] U [-E , 00). 

The proof of Proposition 4.2 is in the appendix of section 17. 
Note that if uo is an eigenvalue of H: 

a 3 Cg =Hg = cug (4.17) 
then uj is an eigenvalue of H* with corresponding eigenfunction a 3 g: 

£cr 3 (a 3 g) = H*(a 3 g) = ua 3 g. (4.18) 

Therefore, we have 
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Proposition 4.3. (1) a(H* ) = a(H ) 
(2) 

W)=span{(!) , (4.19) 

iV(H* T = span {<r 3 ^, ^} = C^} (4.20) 

Here, 

C = °"3^ an d C-m = -^sC-m> ( 4 - 21 ) 

where this choice of £_ M is taken so that (C-^, ^-/u) — 1 111 the l^ol I limit. 

Nondegenerate basis for the discrete subspaces of H and Hq 

For fixed E ^ E , the basis of N 9 (Hq) displayed in (4.19) is a natural basis 
due to its direct connection to the symmetries of NLS. However, this basis is 
degenerate and singular in the limit E — > E *, as we shall now see. 

Consider the basis of N g (H ) displayed in (4.13). Beginning with the first 
element of this basis, explicitly we have: 

-(^)=-(^(:;n'))=^(0 f '- (422) 



where 



and 



a = \a \e™, G = \ " ij0 , (4.23) 



E = ip (-] \a \ 2 ) = Vo* + |«o| 2 x(-, |«o| 2 ) 

= ^o* + |«o| 2 x(-,0) + xi 0) . (4.24) 

Here and subsequently we use the notation x(',p) to denote a generic real-valued 
localized function of x with smooth dependence on a parameter, p, and xi^ is 
localized in x and O(\oij\ k ). 

A nonsingular element of the N g (H ) is obtained by dividing out p . We 
therefore define 

6u = a 3 G ( \ \ F (pl) (4.25) 
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We now turn to the second element of the basis displayed in (4.13). First note 
that 

1^0 = e^-— (K|^o(-;KI 2 )) 



d\ao\ d\ao\ 



V'o + |«o| x(-; Kl ) 

Differentiation of (3.4) with respect to |a | yields: 



(H - E)d laol ^o + 2|* r9|ao|*o + *o3«o|*o = (d M #)*o (4.26) 

Taken together with the complex conjugate of (4.26), this yields, after multipli- 
cation by a 3 : 

a 3 HG ( 1 ) = (d ]ao] E) a 3 G ^ \ ) H^o 

= K|£oi (4.27) 



We define 



£02 = Go ( ! J (4.28) 



where 

d 

F o = oi — T (l Q! o|^b(-; K| 2 )) = ^o* + K| 2 x(^ K| 2 )- (4.29) 

C/|Q!o| 

By the above calculation £02 lies in the null space of (a 3 7i) 2 . Therefore, the pair of 
vectors: £ i and £02 spans N g (H, ) and is nonsingular as E — > E^*. By a previous 
remark: 

C01 = o"3^02 and C02 = ^01 (4.30) 

form a nonsingular basis for Ti^. This choice of basis will facilitate a uniform 
description of the dynamics in a neighborhood of the origin. 

The above construction and Proposition 4.2 imply the following basis for the 
discrete subspaces of Hq and Hq. 

Proposition 4.4. 



N g (H ) 


= span{£ i,fo2} 


(4.31) 


N 9 (K) 


= span{Coi,Co2} 


(4.32) 




= span{^,^_ M } 


(4.33) 


N(H* T lA 


= span{C M , C- M } = -o-s^-m) 


(4.34) 


(C«,6> 


= C ab 5 a b + 0( \a 2 ), 


(4.35) 
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where a and b vary over the set {(01), (02), /x, — //}. For a small we have the 
expansions 




(4.36) 



and = 

Finally, we shall find it useful to note that 

Coi = o"3^02 = foi + \a \ 2 x(x; |«o| 2 )- (4.37) 

Estimates for the linearized evolution operator 
Theorem 4.1. (Linear dynamical stability [54]) Let 

M = NgiHl) 1 - n (H 1 x H 1 ) . (4.38) 
There exists C > such that for any f £ M 

\\e- iHot f\\Hi < C\\f\\ m (4.39) 
Theorem 4.2. (Dispersive estimates) [10] Let 

Mi = [N g (H* )(BN(H'^n)] ± (4.40) 

and P c the associated continuous spectral projection. For any q > 2 there exists 
Ci,g > such that 

\\e- inot P c f\\ L « < C M rt + f||/|| iP , (4.41) 

where p' 1 + q^ 1 = 1. 

Remark 4.1. We shall in later sections use the notation M(t) and M\(t) to 
denote corresponding time-dependent subspaces relative to the time dependent 
operator H^t). 
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Theorem 4.3. (Local decay estimate) Let a be sufficiently large. Let oo G { v e 
IR : \v\ > \Eq\ }, the interior of the continuous spectrum of Ho. Then, for t > 

\\(x)-*e- i ™P c (?i -u-iO)- l (x)-' r \\ B{L 2 ) < C(t)-i, 1 = 0,1 (4.42) 
||(x)-^ fe e-^P c (x)- CT || B( ^ )i2) < C(t)-I~ k , (4.43) 

where H = H + E cr 3 . For t < 0, tie same estimates hold with —iO replaced by 
+i0 in (4.42). 

Theorem 4.3 follows from the technique of [23] and results in [11]; see also [24]. 
In particular, for a Schrodinger operator Ho = —A + V, one has the local decay 
estimate 

\\ {x )-° e -^P c (H )(x)-°\\ Bm <C(t)~l, (4.44) 

where a is positive and sufficiently large. The key to the proof is an analysis of 
the resolvent, (H — A)^ 1 near A = 0. One uses the spectral theorem 

e -itH p( H \ = / e~ itx E'(X)d\ (4.45) 
Jo 

and an explicit expansion of E'(X), which can be expressed in terms of the 
imaginary part of the resolvent. The expansion of the resolvent is valid in the 
space B(L 2 ((x) a dx), L 2 ((x)~ a dx) for a sufficiently large. Time decay is arbitrar- 

3 

ily fast for functions with spectral support away from A = 0, while the t~z decay 
results from the behavior near A = 0. The analogue of (4.43) is an estimate for 
e~ ltH °HQP c (H ), which follows from the expansion of [23] applied to the formula: 

e- itHo H^P c (H ) = r e' itx X k E\\)d\. (4.46) 
Jo 

Finally, this result can be obtained in the matrix case using the approach of [11]. 

5. Decomposition and Modulation Equations 

Consider the nonlinear Schrodinger equation 

id t <p = H<j> + A|0| 2 0- (5.1) 

In the regime of low energy we decompose the solution of NLS in the following 
form: 

0(t) = e- ie « [Mt) + + Mt)] . (5.2) 
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Here, ^o(t) and ^i(t) represent motion along the ground state and excited state 
manifolds of equilibria and 02 is a decaying correction term, lying in an appropriate 
dispersive subspace. The phase, 9 is divided into two parts: 

e(f) = e (t) + e(f), (5.3) 

where 

9 (t) = f E (s) ds. (5.4) 
Jo 

Thus, d t Q (t) = E (t) is the modulated ground state energy, and 0(£) is a "long 
range" logarithmic correction, which is to be derived below 
We begin by setting 

= e- i0 W [¥„(*) + <h] ■ (5-5) 
Substitution of (5.5) into (5.1) yields the following equation for fa: 



id t fa = (H-E (t))fa + 2A|^ o (t)| 2 1 + M&tjfa 

+ (-^e + 2A|0 1 | 2 )^ o (t) + A^(t)0 2 + x\fa\ 2 fa-d t e(t)fa 

idtMt) (5-6) 

Next, we decompose fa into a part along the excited state manifold and a 
correction term: 

fa EE + 2 . (5.7) 

We have, using (5.6), 

^02 = (F- J B o (t))0 2 + 2A|^ o (t)| 2 2 + A^(t)02"-9 t e(t)02 

(^,i(*) + ^9(0+^)^(0 + 2A|^ (t)| 2 ^iW + A^oW 2 ^) 
+ (-d t Q(t) + 2A^ o (t)|0i| 2 ) *o(*) + A^(t)0 2 
+ A(|0 1 | 2 1 -|^ 1 (t)| 2 ^ 1 (t)) 

- (5.8) 

Since equation (5.8) involves 02 it is natural to consider the system governing 2 
and 2 . Let 




(5.9) 



23 



and introduce the matrix linear operator: 



H-E (t)_ + 2\\V (t)\ 2 Att§(t) 

\*l(t) H-Eo(t) + 2\\* (t)\' 



w) - °a 1 owi „ °o;, 1T ,^ l2 • (5.10) 



Then we have (recall that 0i = + 2 ) 

^$ 2 = H (t)<S>2 +F-id t (^ ) - ( (£oi(l) + $e(f))<7 3 + iSt ) ( ^ ) , 

(5.11) 

where 

f ee -$e(f)<7 3 ( *° ) - ^e(t> 3 ( ^ ) 



c.c. 



+ 



5.1. Modulation Equations 

Motivated by the results of Theorem 4.2 on dispersive decay, we shall require that 

$ 2 (f) G A<i(f), (5.13) 

where .Mi is defined in (4.40). Equivalently, 

P c (Ho(t))$ 2 (t) = * 2 (*), (5.14) 

where P c (Tio(t)) denotes the continuous spectral projection of Ho(t). By Propo- 
sition 4.3 this imposes four orthogonality conditions on <3> 2 (t): 

(^a(t),<f 2 (t)) = 0, (5.15) 

where a G {(01), (02), /x, — //}. We impose (5.15) at t — and now derive modu- 
lation equations for the coordinates a (t) and OJi(t) ensuring that (5.15) persists 
for all t 0. 
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To derive the modulations we first take the inner product of (5.11) with the 
adjoint vectors cr^a to obtain the identity: 

(asta , id t ( ll ) ^ + (asta, ( (£01 + d&)a 3 + id t ) ( ^ ) ^ 

= (K(t)^a, $2> + (<7 3 £„, F) + i(d t (a 3 C a ), $ 2 > - i$<<7 3 £a, $2) (5.16) 

The initial data for NLS is decomposed so that (cr 3 £ a (t), $2^)) = for t = 0. 
In order for this condition to persist for all time it is necessary and sufficient that 
the last term in (5.16) vanish, or equivalently: 

Proposition 5.1. The condition that P c (Ho(t))&2(t) = &2(t) is equivalent to the 
following modulation equations for the coordinates «o and a± which specify the 
dynamics along the ground state and excited state manifolds of equilibria: 

, idt^o) + (a 3 C a , (<7 3 (£bi + d t Q) + id t ) *i) , 
' 2|^ | 2 ^i + \\ 



c.c. 



+ /e a ,(-9 t e(t) + 2A|0 1 | 2 ) ( * 



-d t e(t)(S a ,$ 2 )+ i(d t (a 3 Ca),^2). (5.17) 

where a E {(01), (02), /i, -u} and ^ = (^,¥~), j = 0, 1. 

Remark 5.1. (AJ Note that the term {HQ(t)a 3 ^ a: $ 2 ) hi (5.16J vanishes by the 
orthogonality constraint and because H^t) maps the discrete subspace into itself. 
It therefore does not appear in (5.17). 

(B) The last term in (5.17) is present due to the time-dependence of the eigen- 
vectors Ca(t)- An important simplification of this "commutator term", which we 
require for a = (01), is carried out in the appendix of section 18. 

Initial data for the system (5.17), (5.11), governing ceo, oi\ and $2 are obtained 
as follows. Given data 0o for NLS, we find ceo so as to minimize 

||0O-*aol| 2 ; (5.18) 
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see [39]. This a® is used to define the initial Hamiltonian 7Yo(0). Now decompose 
0o using the biorthogonal decomposition associated with Ho(0). This specifies 
a (0),ai(0) and $ 2 (0). 

5.2. Conservation laws and 'a priori bounds 

In this subsection we obtain bounds on ao, a± and 02 using the conservation laws 
of NLS, noted in the introduction. 

By the L 2 conservation law: N[0] = N[0 O ] we have 

h(t) = \* ( t )\ 2 + \* 1 (t)\ 2 + J \Mt,x)\ 2 dx 

= N[0 O ] - 23? | ^ ^7 - 23? J Vofo - 3? / (5.19) 

Note that / = C(|ct | |«i| 2 ) + O(\a \ 2 |«i|). Furthermore, we have 

the orthogonality relation (5.15) with a = (01), (o"3^oi) ^2) = or equivalently 
3? J \I/o02 = 0. Therefore (5.19) is equivalent to 

( 1 + 0(h(t)) ) h(t) = N[0 O ] - 3? I ^i02~ (5.20) 

The orthogonality relation (5.15) with a = \i or (03^, $ 2 ) = implies 3? J \l/i02 = 
O (|ai| 3 + I oil I |«o| 2 )) ||0o||2 and therefore if N[0 O ] is sufficiently small 

|^ (t)| 2 + |^i(t)| 2 + / \Mt,x)\ 2 dx < CN[0 o ] (5.21) 

In a similar way, the conservation law H en [0(t)] = H en [0 o ] implies, if So is suffi- 
ciently small, the bound: 

J |V0 2 (t,x)| 2 dx < C£ (5.22) 

6. Toward a normal form - algebraic reductions and fre- 
quency analysis 

We now embark on a detailed calculation leading to a form of this system, which 
though equivalent, is of a form to which normal form methods can be easily 
applied. 
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6.1. Modulation equations 

Equations (5.17) are a coupled system for a , a± and their complex conjugates. It 
is natural to write the system as one which is nearly diagonal. This can be done 
by taking appropriate linear combinations of the equations in (5.17). 

The equation which essentially determines «o can be found by adding the two 
equations obtained from (5.17) by setting a = (01) and a = (02). This gives: 

MZoi + £02) , id t * ) + (<7 3 (&>i + £02) , (cr 3 (E 01 + d t G) + id t ) ^> 



— A ( £01 + £02, 



2|^ | 2 ^i + ^0^1 



c.c. 



+ ^01 + £02, ( -d t e(t) + 2A|0!| 2 ) ( *° ) ^ 

-9 t e(0(£oi + £02, $2) + i<$(<7 3 (£oi + £02)), $2). (6.1) 

The difference of the a = (01) and a = (02) equations is the complex conjugate of 
the equation (6.1). 

The equation which essentially determines a.\ is equation (5.17) with a = \i: 

(<7 3 ^, idttio) + (<7 3 ^, ((73(^,1 + ^6)+^)*!) 



c.c. 



+ ^,(-W) + 2A|0 1 | 2 )( *° ^ 

-a t 0(t)(^,$ 2 )+ i(9 t ((7 3 ^),$ 2 ). (6.2) 

The equation corresponding to a = —fj, is the complex conjugate of this equation. 
Since 

+ (&)** (6 - 3 » 

we see that the equations for dj = (a , a ) T , j = 0, 1 can be expressed in the 
form: 

iM 00 d t ( ^ J + iMoi^t ( ^ ) + (^01 + <9i0)N O i ( ^ ) = F (6.4) 
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iM 10 d t (21+ 'MiiSi (£ i + $e)Nu f ^ 



Here, for fc = 0, 1: 



N 



01 



M 



lk 



M ok = G 



Nn = 



= F, (6.5) 



o"3(6i + £02), ^= 



^3(^01 + £02), 
03(61 + 62), 




03(61+62), ( |p| 




03(61 + 62), I ^ 



03(61 + 62), 
0-3(61 + 62), 

dk*k \ \ \ 
dk^u )/ 

)), 




'k^k 



6.6) 



(6.7) 



(6. 



(6.9) 



6.2. Algebraic reductions and determination of 0(£) 

To express the modulation equations (5.17) in a tractable form we shall make use 
of a number of notations and relations which now list for convenience; see also the 
section 16. 



xi^ denotes a spatially localized function of order |<x,| fc , as \atj\ — > 0. 



denotes a quantity which is of order \aj\ k as \aj I -> 0. Both xl and O k J ' are 
invariant under the map aij 1— > a^e* 7 . 

oi ^ = 0£>o£\ fc = jb 1 + jb 2 . 



,0) 



0i 



*i + 



a = K|e* 70 , c^e 2i7 ° 



(9 fc * fc ,9 fc *fc) = ^k* + \ak\ 2 xi k \alx { o ] ), k = 0,l 
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U + ^° 2 ~ Go [fo-F>)-[ e -no x (o) ) 

~ { WxP ) 

F + F^ = 2^o* + |«o| 2 x(-;l«o| 2 ) = 2(^o* + X2 0) ) 
F -F^ = |« | 2 X(-; l«o| 2 ) = X^ 0) - 
(^^^(SO)) = 0, (6.10) 

Using (6.10) in (6.1) we get: 

2i(l + O^)d t a + ialO^dtOTv + (<7 3 (£ i + £ 02 ) , ( <r 3 {E 01 + 00) + id t ) ^) 
= -2(^0(F ( ;,M/ o )-2A(F ( ;,|0 1 | 2 M/ o )) 

+ (F + Fq, 2A|^o| 2 ^i + X^Wi + Al^il 2 ^! - A^il 2 ^!) - e 2 ^> (F - Fq, c.c.) 
+*e 2 ^' (d t M&i + 62)], $2) (6.11) 

Determination of 6(t) 

We anticipate that generically |0i| ~ |o;i| ~ for t very large, ao will have 
a limit as t — > ±00 if dtao is integrable. We ensure this by choosing 80 to cancel 
the terms which are of order t _1 and nonoscillatory. Thus, we choose to satisfy 

d t e(F^, * > - 2A(f; i^i 2 ^) = o. (6.12) 

To leading order this gives: 

^0~2A(^,|0 1 | 2 ). (6.13) 

In this way, a logarithmic correction to the standard phase, J* * E (s)ds arises. 
Equations (6.11), (6.12) together with Proposition 18.1 imply 

2i(l + Of ] )d t a + iaftO^dtm + (#oi + 00) (£ i + 62 , *i) + (^3(61 + 62), 
= -^e(xi° ) ) 02)+^0a 2 (x? ) ,0;) 

+A^ 0) (|a | 2 «i + agar) + A(^ 0) , I0i| 2 0i - l*i| 2 *i> - Aa 2 ,^, I0i| 2 0i" - l*i| 2 ^> 



^|«o| 2 (x ( o\<p2) + e 2 ^(x { o\<p2) +Aa \ 2 d t lo (x'o> , fa) + ^ ( X ( o>, 02> (6.14) 



A ) ^70/^(0) 
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We now turn equation (6.2). Using (6.10) equation (6.2) can be written as: 

(1 + Of> l) ) id tai + (O(a 2 ) + 0{a\)) idtW 
+ ( (1 + C^Vi + alaiOfV )(E 01 + d t Q) + (a 3 ^, id t $ ) 
= -d t Q( (xA2)+a 2 (xA2))-d t eol?' 1) a 
-XO^daol 2 ^ + a$ai) + A^ 0) a (x, I0i| 2 > 
+\a^(x,4>\) + A«o(x,0i 2 ) 

+\( X , |0i| 2 0i - |^i| 2 ^i) + \a 2 (x, ~ l^il 2 ^) 

+i(d t (a 3 Q,<£> 2 ) (6.15) 

We next write the systems for c? and 6l\. 

iM 00 d t do + iM 01 d t a 1 + (E 01 + <9i6)o- 3 N iai 
\/n(o) ( l«o| 2 «i + ajjoT \ 



" 1 c.c. 



c.c. 



«M u a t «i + iM 10 d t a + (£ i + ^6)a 3 Nnai 

Z' |ao| 2 ai + «o«7 \ 
- AC?0 a *[ c.c. ) 

+Xa ( (X,|0i| 2 0i-|^i| 2 ^i) + ^(x,|0i| 2 ^-|^i| 2 ^) \ 



c.c. 



+A(J3 ^ «o(x, |0i| 2 ) + MX, <Pl) + «o(x, 0i ) ^ 



c.c. / \ c.c. 



The matrices Mj k , Nj 1; j, k = 0,1 are displayed in equations (6.6-6.9). For 
the (generic) case where we expect | cto j to approach a nonzero limit and a± and 
02 to decay to zero, we shall use the expansion: 
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/ 2(1 + Of ) a 2 of 
M °° " { ^O?) 2(1 + Of) 



1 + Of 1} OK) + 0(a 2 ) 
OK 2 ) + O(or 2 ) 1 + Of 



M 



li 



The matrices Mi and M u are higher order in a and satisfy: 

' of a 2 O, {0) ' 



The matrices N i and Nn satisfy: 



N 



Of + 0f a 2 (Of + 0f^ 
I V(Of + Of) Of + Of 

1 + Of a 2 Of 1} ' 



11 " U^r i+of 



(6.18) 



M 10 ,M 0lH _ 2 ^o) ^ j (6.19) 



Simplifications to equations (6.1-6.2) 

(1) Since Mi and M i are higher order in a we can eliminate d t a\ from (6.16) 
and d t a from (6.17). 

(2) Note also, that "commutator terms" with factors like <9 t | a; 1 2 or d t a.Q can be 
eliminated via redefinition of the near identity matrix M o through incorporation 
of a higher order correction. 

(3) We can eliminate the term proportional to |a | 2 <9t7o as follows. Consider the 
last two terms of (6.16). Since <9 t |o;o| 2 = o:od t ao + a dtao we can incorporate the 
second to last term of (6.16) as a higher order correction to the near identity 
matrix M o, M 00 . Our goal is now to eliminate the eliminate d^o from the 
equation. The system now be written as form: 

id t a = Moo" 1 ( ... + *0(K| 2 )<9 t7 o) (6-20) 
The first component has the form: 

id t a = 1 st component of the vector : 

M 00 _1 (...) +iO(\a \ 2 )d tl Moo~ 1 (6.21) 
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Since «o = |«o|e 170 we have 

i<9 t |aio| — |a |9 t 7o = e~ no x 1 st component of the vector : 
M o _1 (-) +iO(\ao\) |a |<9i7oM o _1 



(6.22) 



By taking the real part of (6.22), for |a | small, we can solve for |a |9 t 7 . This 
enables us to eliminate it from equation (6.16) as a higher order term. 

Implementation of these simplifications leads to the following 
Proposition 6.1. 



id t a 



M 



# 
oo 



-{E i + <9 t e)<7 3 N i<Si 
c.c. 

'■' L,.|2^. _ | lT , |2 lT ,.\ f a 2/, <"l 

C.C. 

(0) A \ , „.2/„(0) 



, \/n(°) / l a o| «i + "o^i 
+AO ° CT3 ^ c.c. 



(xr,0 2 )+^( x ^,0 2 ) 

c.c. 



(6.23) 



idtQ-i 



A(t)c?i 



# 



Acr 3 



( X , |0i| 2 0i - l^i| 2 ^i) + a%(x, |0i| 2 0i - |^i| 2 ^i) 



c.c. 



«o(x, |0i| 2 ) + «o(x, 0?) + «o(X> 0i' 

-<9 t 6 o- 3 



c.c. 



(x,<h) +«o(x, 02) 



c.c. 



(6.24) 



Here, 



&10 + ^0 l a o| ^0 a o 

n^'^Tv- 2 p ro^'^u 1 2 



(6.25) 



A 22 

01 



-A n , A 2 i 

^1 + 02 



-A12, and 



M* and Mfj are near-identity matrices, whose deviations from the identity give 
rise to higher order terms which are subordinate to the leading order behavior 
obtained in the analysis which follows. 
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Finally, equations (6.23) and (6.24) are coupled to the equation for $ 2 given by 
(5.11). 

6.3. Peeling off the rapid oscillations of cki 

Fix T > and large. We rewrite (6.24) centered about the ground state at time 
t = T: 

%d t a x = A(T)a 1 + (A(t) - A(T)) a l + F, (6.26) 

where 



Am / E 01 (T) + O^\T)\a (T)\" 0^\T)o^(T) , 



Since A(T) is a constant coefficient matrix it is a simple matter to obtain the 
fundamental matrix. 

Proposition 6.2. The system 

id t a x = A{T)a 1 (6.28) 
has a fundamental solution matrix: 
X(t) = 



c 21 c 22 /V e 



-i\-t 



cti cr 2 j ^ e- A +* 
1 O^ l \T)a {TfE^{T) \ ( e" lA +* 



O^\T)a (TfE w \T) 1 J V e ^ / 

(6.29) 

The eigenfrequencies, A±(T) are given by A+(T) and A_(T) = — A+(T), wnere: 

A + (T) = ^o^+OrVo^)! 2 , (6.30) 
wiere £io = E\ — E Q , and provided |o;o(^)| 2 /-Eio(^) JS sufficiently small. 
We use the fundamental matrix, X(t), to define a change of variables: 

5i = X(i)0 = X(t) ^ ^ ^ (6.31) 

Therefore, 

ai = e-^^i + O^aK^E^^e-^Ji- (6.32) 
33 



Then, (3 satisfies 

id$ = X _1 (t)F(X(*)/3(*),t) (6.33) 
Note that since the linear in a\ terms have been removed by the change of variables 

(6.3i), m ~ o(iai 2 ). 

Note that by (6.29) 

x = (e^ \ / l + O^\T)\a (T)\^(T) -O^ (T)a$(T)E£(T) 
V -O^ \T)a (Tf l E w \T) 1 + (T)\a (T)\*E w 2 (T) 

(6.34) 

Written out in detail, from (6.33) and the various definitions, we have 
Proposition 6.3. The equation for f3\ has the form 

idtPi = Xe iX+t ( X o, \e~ iX+t Mu + 2 |>o 

+ Ae a +*( Xo , (e" iA+t A^i* + 2 ) 2 )^ 

+ Ae a + t ( X o,(e iA+ 'A^ + 2 ") 2 )«o 

+ A( Xo , [\e~ lX+t Mu + 2 | 2 (e- iA+ *AVi* + <h) ~ e" iA+t |A| 2 /5i<] ) 

+ 7^/3, where (6.35) 

7fy = X- X (t) [(A(t) - A(T))X(t)/3 1 + ^ 1 ] (6.36) 

Remark 6.1. 

X- 1 (t)[A(t)-A(T)]X(t) = Real Symmetric Diagonal S(t) + e - 2a +*B(t) (6.37) 
Therefore, the nonoscillatory part involving S(t) does not effect tie evolution of 

IAI 2 . 

6.4. Expansion of 2 

The next step is to get an appropriate expansion of 2 , which upon substitution 
into (6.35) can be used to isolate the key resonant terms in the f3\ equation. The 
analogous steps are then repeated for the oo equation. Finally, a near-identity 
change of variables is constructed which maps the system for o and fli to a new 
system (a normal form plus corrections) for which the dynamical behavior is more 
transparent. 

$2 solves equation (5.11). We shall require dispersive decay estimates for $ 2 
and these are most naturally obtained relative to a time independent Hamiltonian. 
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Since eeo(t) is expected to tend to a limit as t — > oo and since we are fixing a time 
interval [0, T], it is natural to use as reference Hamiltonian, the operator TCo{T). 
We now make use of the linear spectral theory of section 3 and decompose $2 into 
a part lying in the "discrete subspace of Ho(T)": 

N g (H (T)) © N(H (T) - /i(T)) © jV(ft (T) + fi(T)) (6.38) 

and a part lying in JAi(T), the "dispersive subspace of Tt (T)"] see (4.40). 
Let [3, 11] 

$ 2 = k + n + 77, (6.39) 



where 



fc = ^ (^(T),^ (6.40) 

5aeAf s (H (T)) 

]T (a 3 6(T),$ 2 )6 (6.41) 

fteiV(Wo(T)TM(T)) 

77 = P C (T)$ 2 = $ 2 -A;-n (6.42) 



Since 



= 0, £ a e N g (Ho(t)) 

(<t 3 &(*), ^(0) = 0, ^ e W fl (Wo(*) =F M*)), 

may be replaced by — £(T) in the definitions of and n. Inserting the 
expansion (6.39) into (6.40) and (6.41) and defining 

p Dun (t,T) = ]T (a 3 (£ a (T) — £ a (t)), -)£ a (6.43) 

? eAf 9 (Wo(T)) 

p neut (t,T) = £ -&(*)),•>£« (6.44) 

&eiV(Wo(T)TAt(T)) 



we have that k and n may be expressed in terms of 77 as follows: Then, 



j _ ( Pnuii(t,T) p nu n(t,T) 
{ p Qent (t,T) p ncut {t,T) 



n ) \ p ncut (t,T)r] J ' 



(6.45) 



Therefore, we have 



Proposition 6.4. There exists Eq > such that if \cto(t) — cto(T)\ < Eq then the 
relation (6.45) can be inverted and we have 

$2 = $ 2 W = k[rj] + n[rj] + 77, (6.46) 
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where $2 is linear in 77 and continuous in the weighted (local decay) norm of 
f^\\(x)-°f\\ 2 . 

The statement about continuity in the weighted norm follows from the spatial 
localization of the generalized eigenfunctions. 

Remark 6.2. Note that since |f (T) - £ a (t)\ < C£^\a (T) - a (t)\, we have the 
simple estimate: 

\k\, \n\ < C\a (T) - ao(t)\ ^"^(Olk (6.47) 

Anticipating that fort very large, \ct (T) — ct (t)\ ~ and \\{x)~ a &2(t)\\2 ~ t~^, 
it follows that \k\, \n\ ~ t~ l . Thus, \k\ and \n\, are expected to decay faster than 
77. 

Finally, rj satisfies the following evolution equation obtained from (5.11) by explic- 
itly introducing the reference Hamiltonian, Ti, (T), and applying the projection 
P C (T) to the equation. 



id t T] = H (T)77 

C.C. 

I 2 iTr I lTf 2i 



+ («„(«) - H (Tj) n + (£,(«) - Eo(t)) P c (T)<r 3 



+ AP C (7>3 

+ AP C (7> 3 



c.c. 



1^1 + VTi$ 2 W| 2 (^l + 7Ti$ 2 M) - l^l| 2 *l 



C.C. 



*P C (T) ^ Vl ^ 1 W ' dtSl ^ + V °*°(*) ' 9 '«oW j (6.48) 
7. Normal Form and Master Equations 

In section 4 we decomposed the solution, 0, in terms of coordinates ao(t) and ai(t) 
along manifolds of nonlinear bound states and 2 , a correction which lies in a time- 
dependent subspace, AAi(t), of continuum modes. 2 (£) was then decomposed 
into its discrete (k and n) and continuous (77) components with respect to a time 
independent Hamiltonian, H (T). We also observed that k and n are determined 
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by and are expected to be more rapidly decaying than rj (Proposition 6.4 and 
Remark 6.2). Therefore, the evolution of <j) is determined by a®, a\ and rj. Finally, 
the fast oscillations of a± are removed by the introduction of Pi = [X(t)di]i ~ 
e" iA +'«i. 

We now seek a form of the system for a an d Pi from which the large time 
dynamics can be deduced. We obtain this "normal form" by first solving for rj 
(see (7. 21), (6. 48) ) as a functional of cto, Pi and the initial data r/(0) and then 
substituting an appropriate expansion (see sections 7 and 8) into the equations 
for a and Pi. 

Proposition 7.1. The Normal Form There exists a near identity change of 
variables 

( &o \ = ( ®o\ , ( J a [a ,Pi,t] \ , s 

v Pi ) ~ { Pi ) \ Mao,M J [l > 

where 

J k [a ,p 1 ,t] =O(\a \ 2 + \p 1 \ 2 ), k = a,P (7.2) 
and bounded uniformly in t, and such that 

id t d = (C1022 + ^r w ) |/3i| 4 a + F a [d ,Pi,r],t] 

idA = (cn2i-2ir w )|(5o| 2 |A| 2 A + Fp[a ,Pi,r), *]■ (7.3) 

The properties of F a and Fp are briefly discussed in the remark 7.2 following 
Corollary 7.1 below and described in detail in section 8. 
Furthermore, 

T = r w , +£>(K(T)| 2 ) > 0, where, 

uj* = 2E U -E * (7.4) 

The coefficients Cki mn = Oq*' 1 ^ are real constants multiplying monomials of the 
form aQa l P^ i P 1 . 

Remark 7.1. Due to our choice of the phase correction, Q(t) (see (6.12), a term 
of the form CiquOq 1 ' 1 ^ \Pi\ 2 do, is absent from the differential equation for do in 
(7.3). 
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Now let 

P = |ao| 2 and Pi = \~^\ 2 (7.5) 

denote the (renormalized) ground state and excited state powers. Then, by (7.3) 
we have the Nonlinear Master Equation: 



Corollary 7.1. 



= 2rP Pf + P (7.6) 



= -4TP P? + Pl (7.7) 



dPi „_ n2 



dt 



wnere 



P = P [<5 , A,77,t] = 29f(fi F a ), 

Pi = i2 1 [ao,A,»7,t] = 29f(AF^) (7.8) 

A more precise and revealing variant of Corollary 7.1 is Proposition 11.1, which 
is stated and proved in section 8. 



Remark 7.2. The terms F a and Fp are such that P and Pi are not small per- 
turbations of the leading order terms in (7.6, 7.7) for all t > 0. In fact there are 
three time intervals defined in terms of transition times to and ti (see section 8), in 
which we consider the system (7.6-7.7): I = [0, to], h = [to,ti] and I 2 = [ti,oo). 
It is only for sufficiently large time, (t G h), where Po and Pi are negligible. The 
behavior on short (0 < t < t ) and intermediate (t < t < t±) time scales can be 
very different. We go into the details of P and Pi in section 8 but wish to make 
some remarks at this stage which indicate our approach. 

If we drop the terms Rj then we have a Bow, which evolves in the first quadrant 
of the P — Pi plane according to: 

^ = 2Tp pl (7.9) 

^ = -Wpopi (7.10) 

where solutions for typical data converge to the p axis with a rate (t) -1 . In order 
for the corrections coming from P and Pi to be small, intuitively it is sufficient 
that 

Rj ~^(P Pi 2 + <t>- 3 ), (7.11) 
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where Sq is small and p > 0. This is what we show for t > t\. For the intermediate 
time range, to < t < t\, we show that the behavior is controlled by the system: 

^ = 2rp P 1 2 + Wr 3 ) 

^ = -ATP P^ + O^ Pr) + 0((t)- 3 ), (7.12) 
where m > 3. Therefore, for intermediate times we need to show: 
R ~ ^(PoA 2 + (t)- 3 ) 

Ri ~ ^(PoPi 2 + /^Pr + (i) -3 ). (7.13) 

What makes the analysis subtle is the dependence of Rj on P and Pi in a manner 
which is nonlocal in time. That is, 

V ~ re^ o{t " s) P c (Ho)x«r/3r^ m3 A m4 rfs (7.14) 
Jo 

Local in time terms are simple to dominate by the leading terms. However, 
nonlocal terms require careful analysis. Note in particular, that due to the "history 
dependence" of such terms, being expressed as time integrals from up to t, an 
analysis of the effect of such terms for t > ti requires use of estimates on other 
time regimes t < t and t < t < ti as well. Furthermore, there are no decay 
estimates on either P or Pi in the intermediate interval t < t < t 1 or on the size 
of this interval. 

The normal form of Proposition 7.1 is essentially the Poincare - Dulac normal 
form which can be constructed along the lines explictly implemented in [44]; see 
also [1]. We now give a detailed outline of the procedure with explicit illustrative 
detail of key points concerning the treatment of resonant and nonresonant terms. 



Resonant terms and removal nonresonant terms: Here we illustrate, by 
way of a simple example, how non-resonant terms can be removed by near identity 
changes of variables. Consider the scalar ordinary differential equation 

A'(t) = \A(t)\ 2 e mt (7.15) 

where A(t) is a complex valued function. We shall introduce a change of variables 
A ^ A = A + q 2 (A,A,t), where q 2 (A,A,t) = 0(\A\ 2 ) and q 2 (A,A,t+ ff) = 
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q2(A,A,t), which is therefore approximately the identity for \A\ small, and such 
that 

A>(t) = ^iWpiW + le^liWpI^+^^W,!^)^), (7.16) 

Here, £4 is 2n/fl periodic in t and 

E 4 (A(t),A(t),t) = 0(\A(t)\ A ) (7.17) 

The change of variables can be derived by elementary means. Integration of 
(7.15) gives: 

A(t)-A(0) = f \A(s)\ 2 e ins ds 
Jo 

(7.18) 

Define A(t) = A(i)-(ift) -1 |>l(i)| 2 e <nt . Then, i satisfies the renormalized ODE, in 
which resonant quadratic terms have been removed. The process can be repeated; 
by introducing further changes of variables A 1— > Ai — A+ higher order in A and 
period in t, non-resonant (oscillatory) cubic terms can be removed to obtain: 

A[(t) = l\A 1 (t)\ 2 A 1 (t)+ i k\A 1 (t)\ 4 A{t)+ (7.19) 

where k is real. That the coefficients in the first to terms of this normal form are 
purely imaginary implies that, to this order, the amplitude is independent 

in time. This is the typical situation of the norm form finite dimensional Hamil- 
tonian systems, in which resonances occur between isolated discrete frequencies. 

We next examine resonances between discrete and frequencies and the con- 
tinuum of frequencies, associated with the continuous spectral (dispersive) part 
of TLq. These can introduce nonconservative terms into the normal from (via co- 
efficients with real as well as imaginary parts), which are responsible for energy 
transfer between discrete modes (bound states) and radiation. 

Nonconservative resonant terms and energy transfer: We explain how to 
find the key resonant energy transfer terms, the leading terms in (7.3). These 
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are terms responsible for the exchange of energy among the nonlinear ground and 
excited states mediated by interaction with continuums modes. We focus on the 
Pi equation. Analogous considerations apply to the ao equation. 
Equation (6.35) can be written the following compact form: 

idtPi = E C v^l Jie-™ rt {x Pq r, rf) + 0{rf) + (7.20) 

where C pqr are of order 1, a or higher order in ao, u r E {±A±, ±2A±, 0}, and x vqr 
denote functions which are exponentially localized in space. The equation for a 
has a similar structure. The equation for r\ = r] + r\\ + rj 2 , can be formally solved 
giving: 

V = O(r) ) + 0{r,l) + 0( V 2 ) 

+ E D Piqiri G(^J 1 qi e-^ s x Piqiri ) (7.21) 

Pi,Qi,ri 

Here, G denotes the operator 

/ h-> Gf = -% /"* e -*Wo(*-«) p c f( s ) ds. (7.22) 
Jo 

We insert the expansion for 77, (7.21) into the terms involving inner products 
{Xpqr,v) i n (7.20) and in the corresponding equation for a - This yields a coupled 
system for a and fa which is closed up to higher order. The terms in the resulting 
equations are of the form 

E E C Mr Z) Min ^e-^ ( (v , G^f 91 e^ s X Miri ). (7.23) 

p,q,r pi,qi,n 

We now use the integration by parts lemma: 
Lemma 7.1. 

C e tAs f{s)ds = lim /* e l{A±i5)s f{s)ds 
Jo sio Jo 

= -i(A ± i0)- 1 e iA 7(t) + i(A ± i0) _1 /(0) 

+ i(A±iO))- 1 f e iAs f(s)ds (7.24) 
Jo 

Applying this lemma, we obtain 

e \ Xpqr Xpiqiri / 
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= -e-t^+^ffPi 91 (x Pqr , (H (T) - v ri T ^0)- 1 P c (T)x Pl9l , 1 ) 
+ /3f (0)^(0) ( Xpqr , (H (T) - u ri T iO)- 1 e-^( T ) t P c (T)x Pl , in ) 
+ fl (x Pgr , (H Q (T) - u ri T zO)' 1 e -^o(T)( t -s) e -^ lS p c(T) ^ (ff^^r,) ds^ 

(7.25) 

We first focus on the first term in the expansion (7.25). This first contributes a 
resonant term, which cannot be transformed by a near identity transformation to 
higher order if 

ccv + v ri = 0. (7.26) 

Now consider such a resonant term. We find that in the (5\ equation they are of 
the form: 

-K| 2ai |A| 26l /3i (x , (Ho - u ri =F ity^PcX) (7.27) 

where u ri = —uj r e {±A±, ±2A±, 0}. 

There are two cases to consider: (I) u ri not in the continuous spectrum of Ho 
and (II) v ri in the continuous spectrum of Ho 4 - If v rx is not in the continuous 
spectrum of Ho then the inner produce in (7.27) does not involve a singular limit 
and we get the limit 

(xAHo-Vr^PcX)- (7-28) 

In this case, the coefficient of |cto| 2ctl | 26l /?i is real. Such a term results only in 
a nonlinear distortion of the phase of (3\ and does not effect the amplitude. If u ri 
is in the interior of the continuous spectrum of Ho then the limit is singular. We 
choose the plus sign (+i0) if we study the evolution for t > and the negative 
sign (— iO) for t < 0. This choice related to the condition of outgoing radiation 
explained below; see also [44], for example. Evaluation of this singular limit gives: 

(x , (Ho ~ u ri =F i0) _1 P c x) = A + if, where 

f = 7T ( X , S(Ho - is ri )x) , A = ( X , P.V. (Ho - ^Vx) 

Contributions to the imaginary part are therefore responsible for a change in 

amplitude (here damping of f3\). 

Now, when does a frequency v lie in the continuous spectrum of Ho ? By 

Proposition 4.2, we must have v > —E = \E \ or v < E . By (6.30), A± ~ 

4 In case (II) we consider the generic case where if v ri lies in the interior of the continuous 
spectrum of Ho . 
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±(Ei — Eq). Since v varies over the frequencies 0, ±A±, ±2A± we find that v ri = 
±2A± = —u r resonances are in the continuous spectrum and therefore are those 
giving rise to energy transfer, provided 2E\ — E Q > 0; see (7.4). We now embark 
on the details. 

7.1. Expansion of rj 

We expand rj as follows: 



where r] (t) corresponds to the linear homogeneous evolution with initial data 
r)(0) = P c (T)$ 2 (0) and r/i solves the inhomogeneous linear equation driven by ao, 
ai and i]o(t). 
Equation for rjo(t): 



T)(t) 



Vo{t)+rn(t)+V2(t), 



(7.29) 



id t r] = H (T)r] 
77o(0) = P c (T)$ 2 (0) 



(7.30) 



Thus, 



Vo(t) 



e -*Wo(r)*P c (r)$ 2 (0) 



(7.31) 



Equation for rji(t): 



idtVi = 



H (T) m + P c (T)[Ho(t)-H (T)}$ 2 [r ]o }+ E w (t)P c (T)a 3 




+ 




+ 




2^ |^i + vri$ 2 M| 2 + ^o(^i + 7Ti$ 2 N) 2 



c.c. 





#l+7ri$ 2 [7to]| 2 (tfl + 7ri$ 2 foo]) - |^i| 2 ^ 



C.C. 



(7.32) 



The initial data for r)i is ^i(O) = 0. 
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Remark 7.3. A direct computation using the biorthogonal decomposition of the 
discrete subspaces ofTio and Ti^ of section 4 yields that the second term in (6.48) 
is of a higher order than is explicit: 



E w (t)P c (T)a 3 ^ J = O (K| 2 H + K| 3 ) (xT + X ( o ] + X ( o' l) ) (7-33) 

for | ct 1 << 1 and |«i| I 0. 
Equation for r) 2 (t): 



id t V2 = n (T) m + P c (T)[Ho(t)-n (T)}^ 2 [r 1l + r l2 } 

+ AP (T)a 3 ^ l^ 1 + 7r i $ 2N)vri$ 2 [?7i + T] 2 ] + cc] ^ 

+ AP (7> 3 ^ 2 *oki^ 2 [«i + rf + 2^(#! + Tr^sM)^^^! + n 2 ] ^ 

+ AP c (TV 3 ( ^ (7ri ^ 1+ ) 

+ AP c (T)a 3 2 l^i + 7r i $ 2M| 2 7r 2 <fhi + ^] ^ 

+ AP(T)<7 3 ^ ( 7Fl$2 ^ 1 + ^(^l + TrT^M) + TTl^^l +^ 2 ](^l +VTi$ 2 M) 2 ^ 

+ AP (T)a 3 2 ^ 1 + ^^[^Dki^^i + V2] I 2 + ki^2[«i + rfvri^i + ^2] ^ 

(7.34) 

We expect that n 2 = 0({t)- 1 ). Let 

V2 = V2a + V2b- (7.35) 
By construction we will show that r\ 2a = 0({t)~ l ) and r] 2b = 0((t)~^). 

idtV2a = n (T) V2a + P c {T)[ H (t) - HojT)]^] 

+ 2AP C (2> 3 ( ^o(^ivri$ 2 [^i] +^7r 2 $ 2 [^ 1 ]) + ^ok$ 2 M| 2 \ 



cc. 



+ AP (7> 3 ^ 2 ^07Tl$2hl](^l + 7Ti$ 2 [?7o]) + ^0 (^$2 M )' 
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+ XP (T)a 3 y 2 '^ 1 + ^MIV^M + (*i + 7ri$ 2 M) 2 7ri$ 2 M j 

+ \p c (T)a 3 ^ 2 (*i + 7r i $ 2M) ki$ 2 MI 2 + (*7T^M) K^M) 2 ^ 

+ AP C (7> 3 ( l^ 2 ^ $2 N ) (7.36) 

7.2. Normal form and master equations 

Using (7.33) and explicitly inserting in (7.32) the representation 

ttl = (X(t)/?(t)) 1 = e-^( T ^ 1 + ^ ' 1) ^(T)^ ro 1 (T) e - A -( T )^ (7.37) 
gives the following equation for r\\. 

idfm = H (T) Vl + P c (T)(Ho(t)-H (T))<5> 2 [ Vo ] 

+ E 10 (t)P c (T)a 3 ( |« | 2 Xo 0) + 0i\ 2 X? ) 

e -A + (T) <A + o(°' 1 ) a 2 (T)jBro i (T)e -a_^ \ 



X 



+ AP c (T)a 3 <^i 



c.c. y 

2|a | 2 /3 1 e- iA +( T ) t + a^A + (T)t 



c.c. 



with initial condition ^i(O) = 0. 

Substitution of r/i into the equations (6.23) and (6.35) for ao and Pi gives rise 
to terms which make explicit the resonant exchange of energy between the ground 
state and excited state. We next isolate the key terms in the expansion of rji 
relating to this energy exchange. 

Let's begin with the /3i equation, (6.35). Written out in greater detail we have: 

idtfc = 2A^o*,<)|/5i| 2 «oe lA+{T)t + 2A(^ 1 2 *^2$ 2 )A«o 

+ H^MPfe^™ + 2A<^oA^i*2>C3ia e 2iA +W t + /») 

+ 2X(Mu, |vr 1 $ 2 | 2 )a e iA+(T)t + X(Mu, (vr 1 $ 2 ) 2 )o^e a+(T)i 

+ A(^,7Ti^2)|/3i| 2 e iA+(T) * + A«, {^ 2 f)Jie 2lX+{T)t 

+ A(^,vr 2 $ 2 )/? 1 2 c-^^ + 2A(^,|7r 1 $ 2 | 2 )Ae-^( T ) i 

+ A(V>i*, \^ 2 \ 2 ^2)e lX+iT)t + Tip, (7.39) 
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where 1Z@ is defined in (6.36). 

We claim that the key term in (7.39) responsible for energy transfer is the 
term 

7r 1 $ 2 )A;aoe 2iA+(T)t (7.40) 

on the second line of (7.39). To see this we decompose rji into "resonant" and 
"nonresonant" parts 

Vi = Vir + Vinr (7.41) 

The part rjiNR gives rise to the key resonant energy transfer term in the equation 
for Pi and is the solution to the initial value problem 



idtViR = H {T)i]i R + XP c (T)a 3 ifj^ip 0if aoPie~ , 



j2 e -2i\+(T)t ^ J 

WO) = 0. (7.42) 
Solving (7.42) using DuHamel's principle we have 

WO = -iX /V^W^ftiTl^L^A^le-^t^ ( J ) <fe 

(7.43) 

Recall that the continuous spectrum of Hq(T) is given by points oo such that 
M > l-^ol- By Proposition 6.2, for |a | 2 /-E'io sufficiently small 

\ + (T) = Ei (T) + oS , ' 1) \a (T)\ 2 . (7.44) 

By the hypothesis (7.4), 2Ei* — _E * > 0, if | cko J is sufficiently small then 2A + (T) 
lies in the continuous spectrum of H,o(T). Therefore, (7.42) is a resonantly forced 
system. We expand the solution as follows. Let 5 > and set 

WO = -iAe^^jre^m-^^'PcCTO^^^^W ( J ) **■ 

(7.45) 

Then, r] m = lim 5 ^ Vm- 

We now apply Lemma 7.1 to r]i R with A = Hq(T) — 2A + (T). The result is 
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Proposition 7.2. The limit lim^o Vir — ViR exists in S' and 
ViR(t) = -Xe-^^Mtj^t) (H (T) - 2A + (T) - ^O)" 1 P C (T) ( J ) <Vo* 

+ A^(0)/3 1 2 (0)e-^( T ) t (Ho(T) - 2A + (T) - ^O)" 1 P C (T) ( J ) <Vo* 
+ Ae-*" ^* (Wo(T) - 2A + (T) - ^O)" 1 • 

= + ?7ii» + ^iij c (7.46) 

We substitute (7.46) into the key term (7.40) in (7.39). 

2X(Ml, 7r 1 $ 2 M)A«oe 2iA+(T) * = 2X(Ml, 7r lVlRa )J^oe 2lX+{T)t + R. (7.47) 

Here, R denotes rapidly dispersively decaying terms plus higher order terms in 
|a /9i|. By Proposition 7.2 and the Plemelj formula (16.3) 

= -2\ 2 (^ ^i,n l (H (T)-2\ + (T)-tO)- 1 P c (T) ^ J )^ ^ |<*o| WA 
= -2(A + ir)|a | 2 |A| 2 /5i, (7-48) 
where (using that E - 2E 10 = E - 2(E 1 - E ) = + O(\a (T)\ 2 )) we have 



A = C>J 0) 



= A 2 (^ *<,P.V.(/f-^)-Vo*0 Of (7.49) 

r = r^ef 

= \ 2 n(MlAH-u,)Ml)0^ (7.50) 

Recall that denotes a term of the form l + (9(|a | 2 )- Returning to the equation 
for fa we have: 



id t fa = (A^-ir^)\a \ 2 \Pi\ 2 Pi + .... (7.51) 

We now seek the key terms in the a equation, (6.23). Using that 4>i — + 02 
and the representation (7.37), we have 



2id t a = A(t<)e"^ 



a 
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+ +2A(^L,02)|/3i| 2 

+ A^i*, <PlWie tX+{T)t + \(i>o*i>l,fc)fte- 2iX+{T)t 

+ 2A(^ ^i*, \<p2\ 2 )Pie- a+{T)t + A(^ *, |0 2 | 2 2 ) 

+ K Q0 (7.52) 

We first focus on the key resonant term in (7.52) which is responsible for the 
system settling onto the nonlinear ground state. We claim this term is: 

A^o^L^iV 2 ^* (7.53) 

In analogy with the previous calculation, we have 

A(^o^L,^)/9 1 2 e- 2lA+(T)t = A(^L^2^i Ra )/3 1 2 e- 2a+(T)t + R ) (7.54) 

where R is as above. Therefore, applying Proposition 7.2 we get 

\{^l,^2)Ple- 2lX+{T)t = {~K>* + i^u>J\Pi\ 4 oc + R (7.55) 

The expression R denotes terms which are higher order, oscillatory type and 
dispersively decaying. See also Remark 7.2. 

In summary, we have the following system for a and (3±: 

Proposition 7.3. 

id t a = +«r w J| / 3 1 | 4 a + R 

idtfc = 2(A w .-ir w J|a | 2 |A| 2 A + R (7-56) 

The proof of Proposition 7.1 follows by constructing an appropriate near- 
identity change of variables, transforming (7.56) to (7.3). This is implemented 
as in [44]. 



8. Stability analysis on different time scales - overview 

In Corollary 7.1 we obtained coupled power equations or nonlinear master equa- 
tions governing the (renormalized) ground state and excited state square ampli- 
tudes: 

P = |«o| 2 and Pi = |/?i| 2 (8.1) 
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If we neglect the correction terms Rq and Ri, in (7.6-7.7) we obtain the simpler 
autonomous system of differential equations: 

^ = 2r PoP ; (8.2) 

f - -VM (8.3) 

Note that this system is exactly solvable. Addition of twice (8.2) to (8.3) yields 
that along any solution trajectory: 

2p (t)+ Pl (t) = 2p (0)+pi(0). (8.4) 

This relation can be used to eliminate p\ from (8.2) or p from (8.3). po(t) and 
Pi(t) are thus obtained by quadrature. The dynamics of this finite dimensional 
reduced system anticipates that an initial state, arbitrarily close to but not exactly 
on the excited state branch, with energy distributed among the ground state and 
excited state, will evolve to a state with an increased ground state energy and no 
energy in the excited state. While not strictly correct, since there are nongeneric 
data giving rise to solutions which converge to the excited state [50], this captures 
the generic very large time dynamics. The correction terms R and R\ in (7.6- 
7.7) lead to different transient behavior which may be quite different from that 
suggested by the system (8.2-8.3). However, we show that eventually (t > t x ), this 
system dominates. Moreover, a large class of data, for which the system (8.2-8.3) 
controls the behavior is that for which -Po(O) > -Pi(O) and sufficiently small initial 
dispersive part. 

Before embarking on the details we give a brief overview of the strategy. Using 
the change of variables (a, (3) h- > (a,/3) of Proposition 7.1 we have transformed 
away all local in time nonresonant terms. This introduces contributions to F a and 
Fp, and therefore contributions to R and R\ in equations (7.6-7.7), which are of 
two types: 

(i) local in time terms depending on ao and $i, which can be absorbed by the lead- 
ing terms in (7.6-7.7), with a small correction to the coefficient to the coefficient 
T and are of order b (t)~ 2 ; see Proposition 9.1 below. 

nonlocal in time functions of a and p 1 defined in terms of i] = r][r]Q, rji, 772] in F a 
and Fp. These contribute terms to (7.3) with the same (anticipated) time-decay 
rate as the leading order terms in (7.3). Correspondingly, there are nonlocal in 
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time functions of a® and f3\ which contribute to Rj in equations (7.6-7.7) which 
are of the same (anticipated) decay rate as the leading order terms in (7.6-7.7). 
The goal is to control these nonlocal terms, to the extent possible, by the leading 
order terms. However, due to the different behavior of do and f3\ on different time 
scales the argument is somewhat tricky and we now explain our strategy. 
Let 

to + 1 = snp \r : < r' < r, P (r') < ^ m]x + % ) (8.5) 
r>o { 2 (r - 1){t'Y) 

Propositions 11.2-12.1 and 13.1 will justify this choice, by implying the inequalities 

(8.7). If t < oo, then we have the bound: 

Consider the system for Po(t) and Pi(t), (7.6-7.7). Decomposing Rj into local 
and nonlocal in time parts we have: 

Rj = R l j {P ,Pi,m,r ] ) J r f o K{t,s)rf{s)ds 
= R^Po, P u 770, V) + fj Kit, s)rf{s)ds 
+ /* K(t,s)r?(s)ds 

Jto 



We have that 



4^0,^,770,77) < 0(£S°) ( PoPi + 

^(^0,^,770,77) < 0{E£)( P P? + -^ + O[JP P™) ), (8.7) 



(t) 2 

where £ is the initial total energy and pj > 0. Therefore, for £ small 

w - ^p p?-^+iy (t,sK i ( S )d S 

dPl < -4r'P P? + ^ + j^K 1 {t,s)r^\s)ds + o(jp Q Pr) 



dt ~ u 1 (t) 2 Jt 

with T' ~ T. The reverse inequalities hold with V replaced by T" = T + o(T). 
Next (Proposition 9.2) we introduce the auxiliary quantities 

Q = P - koit)' 1 and Q 1 ^P 1 + h(t)-\ (8.8) 
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where ko(bo, bi,£o) and ki(bo, b±, So) are chosen appropriately and derive equations 
of the form 

^ > 2T"Q Q 2 1 + J t y o (t, S y o d (s)ds 

^ < -4r'VoQl + l t K 1 (t, S )r?\ S )ds + o(jQ~oQ?), (8.9) 

where r" ~ T' ~ T. 

We then proceed with the following continuity argument. At t — t , 

dQo(t ) „ 2 dQi(to) 



dt 



> 2T"Q (t )Qi(t ), < -AT"Q Q {U)Q\{t Q )+0 (JqMQ? {to)) ■ 

(8.10) 

Therefore, by continuity, the following inequalities hold for some time interval 
to < t < t ,i, with r" replaced by 

dQ ° > T'lQoQl ^ < ~2T"Q Ql + O (Jq~oQ?) . (8.11) 



dt 

Let t* = sup{t > t : inequalities (8.11) hold }. We show that 

j t y i (t,s)r n l l (s)ds < 0(£g) ( QoQl + ^ ) (8.12) 

and therefore, up to renormalization of Qj (adding higher order terms of order 
So°kj{t)~ l to the definition of Qj))- Use of this estimate in (8.9) implies (8.11), for 
Sq sufficiently small. The argument can be repeated and therefore, t* = T. 

9. Finite dimensional reduction and its analysis on different 
time scales 

We now begin our study of the generic case, where t < oo and the solution 
converges to the nonlinear ground state family as t — > oo. The following three 
propositions concern the various time-scales which enter the analysis. The first is 
a basic result, a normal form, which is the point of departure for our analysis on 
all time scales. 

Proposition 9.1. Let m > 4. Let 

bo = (*o} -1 (Mx + c^ 2 ), (9.1) 
h = (to)' 1 * ( MI + d£l ) , (9.2) 
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for some order one constants c* and d* . 
If for some to, positive and finite, 



Po(to) > 3 '° (t ;; [ '° ]x) ! (9.3) 



then for t >t 



^ > 2(1 - WPoA 2 - |_ + j Q (9 .4) 
^ < -4(1 - 5 1 )FP P^ + O(Jp P?) + ^ + J U (9.5) 
where Jo and J\ are nonlocal in time terms, which have the form: 

ct 

'to 

The terms encompassed in Jj are derived and estimated in the coming sections. 



J, = f K(t,s) rf(s) ds. (9.6) 

J to 



Remark 9.1. The reverse inequalities of (9.4- 9.5) hold as well with a different 
constant S 2 ~ Si. 

The proof of Proposition 9.1 will be given following the estimates on the re- 
mainder terms, Ri (Proposition 11.1). 

Proposition 9.2. Assume that t is positive and finite as in Proposition 9.1. 
Then, there exist ko = ko(bo,bi,£o) and k\ = ki(bo,bi,£o), such that for t > t 
the auxiliary functions 

Qo(t) ee Po(t)-^ y Qi(t) ee Pi(0 + |y (9-7) 



satisfy 



dQo > 2T'Q Ql + J ■ Elc * 



dt ~ u (t ) (t) 2 



1 < -WQoQi + O(jQoQT) + Ji - ; 7! ' ■ ( 9 -8) 



dt ~ — {to) t {t)2 

for some positive free constants c* and d*. In particular, Qo(t) is monotonically 
increasing for t >t . 
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The next result shows that c* and can be chosen to control the terms Jj. 

Proposition 9.3. (Monotonicity of Qq for r >t ) There exist c* and of order 
one, such that for t > t . 

dQo 



dt 
dQ 



dt 



1 < -4T'Q QI + O(Jq' QT) (9.9) 



The above three propositions are established in the next two sections. We com- 
plete the current section by working out the consequences of the finite dimensional 
reduction (9.9). 

The next proposition, shows that even if Q very small at some stage, it will 
eventually become large relative to Q 1 and the 0{^/Q^Q™) term in (9.8) will 
become negligible. 



Proposition 9.4. Assume t >t and suppose for some r > that 

Q 
Qi 



7T ^ £ o- ( 9 - 10 ) 



Then, 

7^7-7 is increasing for t > t (9.11) 
Qi{t) 

and there exists t x such that 



Qo(t 



Qi(ti) 



4 = ?o- (9-12) 



Furthermore, for t > t\. 



dt 
dQ 



dt 



> 2T'Q Qi 

1 < -IT'QoQl (9.13) 



Finally, for t >t±: 



Ql(t) ~ 1 +4T"g 1 (t 1 )(S [ ^ IQoI) • (* - ti) (9 ' 14) 
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Proof of Proposition 9.2: Define 

Q = P -k (t)-\ Q 1 = P 1 + k 1 (t)-\ (9.15) 

where ho, k\ > are to be appropriately chosen. Note that 

Qo(t) < P (t) and Q^t) > Pi{t). (9.16) 

Using (9.4-9.5) and some estimation we deduce a simplified system for Q and Q±. 
We calculate, omitting the terms J and J 1; which are carried along passively. 
We begin with Q . 

dQ dP k 

~df ~ ~dT w 

> 2YP Q P 2 -^+ k ° 



2 



(t) 2 (t) 



(t)J V (t)J (t) 2 (t) 

> 2TQ Q 2 - AT^-Q Q 1 + 2r A ' : 



2 



2 



We estimate the second term on the right hand side as follows. For any s > 0, 

-4r^g Qi > -2sTQ Qj - 2rg ^. (9.18) 



Therefore, 

dQ 



> 2T(1 - s)Q Q\ + (2TklQ Q (l - -^j + k - fe ) (9.19) 



dt - v ' v x ^ u v s ; u v 

fci > h. (9.20) 



Now set s = t±j and assume 



Then, using that k\ = 10&i, we have 

^ > 2 • ^rg g? + (k -b - isr^g ) (9.21) 

if 

A; = 6 + i8r^ sup g + -2^, (9.22) 

t <t<T \to/ 
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we have by (9.21) and (9.28) 

d -w^-^ + ww (9 - 23) 

In particular, Q is increasing for t >t . 

We now turn to Q\, we have 
dQi dP\ k\ 

= - 4 K g ° + l)( gi -l) 2+ ^ + ^ pn 

= -ArQ Ql + 8r^Q Q 1 -ArQ ^ + b -^^ + o(^p pn 

(9.24) 

The second term on the right hand side is estimated as follows. For any r > we 
have, since 2ab < ra 2 + r~ x b 2 , that 

8r | 0oft <™ + !™i. (9 , 5) 

(t) (t) 2 r 

Therefore, 

^ < -4r(l - ^)Q Ql + 4r^Q ^ + + O(y/P l?) (9.26) 

Let 

E 2 d 

r = 20A;iand6i = ki(l - IQTkiQo) ^4, (9.27) 

(*o) 5 

which is consistent with the constraint (9.20). Then, 

'it £ - 4 'I r ^-JV + 0( ^ Pr) ' (9 ' 28) 

By (9.8), Q is increasing for t > t Q . Therefore, 

1 k 



Qo(t ) > 



100 (t ) 



Po(t) < Qo(t) + ^- = Q (t) + 100 / '° 



(t) y 100(t) 

< Q (t) + 100g (to) < IOIQoW 
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Also, by definition P\ < Qi and so 



dQi 
dt 



< 



19 / — 
-4 • ^QoQl + O(^Q QT) 



(9.29) 



This completes the proof of Proposition 9.2. 
Proof of Proposition 9.4: 



Jt (f) = °° Ql ~f lQ ° ^ ^ { r <?oQ?Qi + 4r ^? - CQTQr} 

> -L { 2 rg [q? - cg^g^gr 1 / 2 ] + 4rg*g 2 } 

Vl 

^ < £q then implies that 



d_ /Qo 

dt IQi 



> 



{2rg [g?-^ /2 gr +1/2 ]+ 4r ^} 



> r 



for m + | > 3 and (g^ < £ << 1. Hence 



(9.30) 



go 

Qi 



> 



Qo 
Qi 



exp(r f Q\{s) ds). 
to J to 



Since Q\ > 0, either Qi { or ff- grows exponentially with t > t . In either case, 



there exists ti > t , such that for t — ii, ^ 

Qo 



ti 



£q. Now whenever, 



— Kb 



(9.31) 



we have 



dQi 
dt 



< 

< 
< 
< 



-AFQoQl + +o( v /g g-) 
-4rg g 2 + Q Qr 1/2 £ r/2 

ATQo [Qj - CS Q - r/2 S^ /2 Q\ 

-4rg g 2 



(9.32) 



56 



for m > | + |. Therefore, by (9.31) we have ^ < 0. Since Q is increasing for 

*i 

t >ti > t , the inequality ^ > £q persists and (9.32) holds for all t > t x . Hence 

< -WQoQl and ^ > 2T'QoQl (9.33) 
or at 

Finally, for t > ti 
Solving this scalar inequality 

Qlit) ~ i + 4r''g 1 (S| ( g ) (t 1 )|(t-t 1 ) (9 ' 35) 

Since Pi < Qi, -Pi(t) decays to zero like (t — h) -1 . This completes the proof of 
Proposition 9.4. 

10. Decomposition and estimation of the dispersion 

In this section, we revisit the decomposition of the dispersive part, 77, which satis- 
fies equation (6.48). Here, we decompose 77 in a manner suitable for consideration 
of the solution on the various time scales. 

Proposition 10.1. 

V (t) = no(t) + eo(t) + Vb(t). (10.1) 

The three terms can be described as follows: 

(i) ?7o(0 is a dispersive wave generated by the data 77(0) = rj m : 

id t T] = Ho 770, 77(0) = Vm, (10.2) 

(ii) eo(t) is driven principally by rjo(t): 

id t e = H e + P c S^[e ,r lb ;r lo \, e (0) = 0, (10.3) 

and 

(in) r) b (t), which is driven by bound state dynamics: 

id t n b = Ho rib +P c S (b) [e ,v b ;m], ^(0) = 0. (10.4) 
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We display expressions for and with the detail required in our analysis. 
We let x denote a generic exponentially localized function of position, x. 

S<°> ee S^ + S[% + S?el + el 

~ (al(t) - «q(T)) xVo + vl + a aiXVo + ®oXvl 

+ (a «i + «?) Xeo + («o + oil) x{Vo + Vb)e + VoVbeo 

+ (<4(t) - «o( T )) X e o + («o + «i)xeo + (?7o + Vb)el 

+ el (10.5) 

~ - ag(T)) 

+ *P C (T) ^ ^i^ 1 ^) ' + Vo*o(f) • ^cfo(t) ^ 

+ (Mt) - E (t)) P c (T)a 3 ) 

+ ( |a | 2 |«i| + |a | |ai| 2 ) x 

+ a a 1 x r r]b + (a + an) xVb + (Vo + e ) nj; 

+ \ao\xVoVb + (vo + e o)Vb + Vb ( 10 - 6 ) 

By (5.22) we have the following H 1 bounds on e and i] b in terms of one 
another: 

||e (*)||tfi < £o+\\vo\\m + \\Vb(t)\\m (10.7) 

\\Vb(t)\\m < £o+\\m\\m + \\eo(t)\\Hi. (10.8) 

Using DuHamel's formula, both the eo and r/b equations can be written as 
equivalent integral equations: 

e (t) = -i f e - iHo ^ P c S {0 \s) ds (10.9) 

JO 

n b (t) = -i f e- lHa{t - s) P c S {h \s) ds (10.10) 

JO 

Therefore, in both cases we must estimate an expression of the form: 

w(t) = f e - lHQ ^ P c S{s) ds. (10.11) 
J 

For estimation, we shall require the following class of dispersive estimates: 
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Proposition 10.2. (1 ) Let 2 < p < oo, p',q > 2,q' and s be related by: 

3 (- - -) = s, p' 1 + (pT 1 = 1, q' 1 + (qT 1 = 1- (10-12) 

Then, 

\\e- iHot P c f\\ P < |*r 3( '-« ) <*>- S( '-' ) (||/||p' + \\d s f\\ q >) (10.13) 

C2) Assume g > 2 and s > 3/q. Then, 

lle-^Pc/IU < |i|- 3 ^-^><t>-f ( H/IK + ||97|| g (10-14) 

Proof of Proposition 10.2: We use the classical Sobolev inequality for functions 
denned on 1R 3 : 

||/Hp<C||07|| 9 , 3(q~ 1 -p- 1 ) = s (10.15) 
and the L p - L p ' estimate (Theorem 4.2) 

\\e' mot PJ\\ p < C hp rS+i ll/H,,, p' 1 + (p'y 1 = I- (10.16) 

For |*| > 1, we use (10.16). For |f| < 1, we have by (10.15) that 

\\e- iHot P c f\\ p < \\d s e- m ^P c f\\ q 

< c \\nle- inot p c f\\ q 

< C \t\~^-~^ \\D s f\\ ql , (10.17) 
q- 1 + (qT 1 = 1 

From estimates (10.16) and (10.17) we obtain (10.13). Estimate (10.14) is ob- 
tained similarly. This completes the proof of Proposition 10.2). 

We now apply Proposition 10.2 with q = 4 and s = 1 > 3/4 to the integral 
equation (10.11) and obtain the bound: 

IMOIloo < f - 1 ,, . 1 ~s {\\S(s)\\ 1 + \\dS(s)U ) ds (10.18) 

We shall use this bound as the first step in estimating ||e (t)||oo and ||?7&(£)||oo- 

More specifically, our estimation strategy seeks a closed system of inequalities for 
the following 
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norms of e : 1. M^j, 2. [de G ]q^.z, 3. [e ] H i. 

norms of r) b : 4. [?7f>]oo;0, 5. [r] b }m;0, 

in terms of norms of the initial data [0(O)]x- 

1. Estimation of ||e ||oo 

By (10.18), we have 

lleoWIU < f - 1 pr t~ 1 X 3 (ll^ (0) (5)lli + ll^ (0) (5)ll|) (10-19) 

so it sufficies to bound and ||^°)|| 4. 

\\S^%: For any t > 0, 

< C(s)-I ( £bfa>(0)U + >/5fo>(0)ft ) < C(£ o )D(rj(0)) 

< C(s)-5 [eoj^.3 ( £ + ^Noo;0 + ^Noo;0 + M^f ) 

+ C(s)-i [eoL.ahoL^WooiO < C(s)-i 

< C(*r 3 [e ]L ;f ( CV^+ Moo;| + Moo;0 ) 

< C(s)- f [eo]oo ; |[eo]^ ; o (10-20) 
|| ftS^ || 4: For any t > 0, 

< C(S)-I (fb[9l»]oo;i + 4Moo;fM^;0 + 4 N^a ) 

< ( ^oMoo;f + ^oNLf oc ;|) 

+ ( £o 5 Moc;oMoo;§ + ^0 Noo.f [l^iO + E Q Moo^Wi^.S^ 

< C(a>-i [eoL ;f ( 4 Noo ; f + 4 M^o) 

+ C(s)-3[e ] oo .| (Mffl;0 + Mh1;o) Mf^O 
3 I 

+ ML; 3 M^jO (M^;0 + M//1;0) 

< [eo]oo ; | [e ]|i ;0 (10-21) 
Use of the bounds (10.20) and (10.21) in (10.19) implies 



USflli 
ll^ 0) eo||i 

ll^lli 
lleolli 



||^ 0) e ||| 
ll^||| 

WdeU! 
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Proposition 10.3. 

IMOIloo < (r f (c(^,ho(o)]x) + £lh(t)Uo + [m(o)U MOW [e L ; | 

+ (*)"§ ( £* + [ m (0)] X + [e (t)] 2 Hh0 + M*)L;f ) ( Noo;f + NL;| ) 

+ (t)' 1 ( Moo;f [eo]^;o + ^o[eo] LL; | + 4[e ]L ; f ) (10.22) 
2. Estimation of [de (t)] L 2 .3 



|<9e (t)|| L 2 oc ~ \\xde {t)\\ 2 (x localized) 

\\x(n )h (t)\\ 2 ( 9 (Ho)' 1 * e B(L 2 ) 

f \\x(n )h- in ^P c S^(s)\\ 2 ds ( by (10.9) ) 

JO 



ft 

< 



(10.23) 

For the purpose of continuing the estimation, we regard as consisting of terms 
of two types: (i) terms having spatially localized (exponentially) functions of x 
as a factor, coming from tyj, j — 0, 1, and (ii) terms like e[j or 7/67/oeo and others 
which are not of this type. It is convenient to refer to such terms as S^q C and 
Snloc below - From (10.23) we have: 

\\de (t)\\L L < J^\\x(n }h- m ^P c S^ c (s)\\ 2 ds 

+ (C + H) Ilx(^o)^-^^)p c 4°i!oc(s)ll2 ds 
* C 



< I - IK*)'SJM ds 



(0) 

(t-s)2 

t-1 ,„ ft 



+ l~ ll e ^ 0( ^ )p c W^ffioc(5)lloorf s + / \\(n )h-^-s) PcS ^ {s)hds 

Jo Jt-1 



Jo (t-s)* 



+ 



(t-8)i 

C r 1 -^—r \\ d S§loc(s)\\oo ds + f \\dS§l oc (s)\\ 2 ds 

JO \t-s\2 Jt-1 

A + B + C (10.24) 



We estimate the integrands of B and C; that of A is similar. 



61 



integrand of B | < C (t) § ( [e^s + [^.3 ) ( [e ] 2 H1 . + [r] b } 2 m . Q + [»)o]|i ;0 

(10.25) 



integrand of C | < (*)"§ ( [e ]^.| ( M/z^o + [%]fli ; o ) 

+ Mooif Moo;|Mffi;0 

+ Mw^~;§ (Moc;0[e ]//i;0 + MoojfM^iO + Noo;|N^;o) 

+ M^i,=c.a[eo]^ ; o ) 



Substitution of these bounds into (10.24), we have: 
Proposition 10.4. 

\\de {t)\\ LlA < C(t)~i ( [e L ;f + [ V o(0)]x ) ( bUf + M^ ; o ) (10.26) 

3. Estimation of [e (t)} H i ;0 

\\eo(t)\\m ~ II^MOIh < [' \\dS^(s)\\ 2 ds (10.27) 

Jo 

Wds^h < c( s )-i (^0N^-;|+4NL ; | 

+ [Vo\m;0 NL;f + £ M W l,oc.| ) 

ll^eolb < C( S )-i ... 

||^ 2) e*|| 2 < C(s)-I [sheo]l,s+sheoL,i[m]m-,o) 

+ NL;f (Mifl;0 + [%]ffi;o) 

+ C(s>"5 Moo;!^]^!,-.! ([r] ] 2;0 + [ribhfi) 
\\de 3 h < C(s)- 3 [e ]L ; | M^o 

Proposition 10.5. For t > 0, 

||e o (0ll^ < C(£b, foo(O)U) + [eoL ; s ( M0)U + M^ ; o ) (10.28) 
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We now turn to the estimation of ||^7&(t) ||oo and ||?7&(t)||,f/i. By (10.18) we have 

MOIloo < f - 1 x a + ds, (10.29) 

Jo \t - s|i (t - s)4 v 3 y 

where 5^ is given by (10.6). 

Proposition 10.6. In terms of Qj, j = 0, 1 we have 

S<« = S^ + sPn b + S?rj 2 b +rj! (10.30) 



where 



# ~ QlQix + QoQlx 



fh\ — — — 

Si ~ Qo Qi XVb + (Vo + e oH + Qo XVoVb 
S?V 2 b ~ (g| + Q?) X776 + fab + eo) r)i 

We now proceed with estimates of in L 1 and in W 1 '^. 
Beginning with HS^Hi we have: 

\\s(% 

< c(q|Qi + QoQ? + QlQiWVbWoo + lklloo(||%||2 + l|eo||l) 



+ cj^holUMh + (Oo 5 + QFJIM« + W»( IMI \\mh + IMh IM| 2 
+ CIMUMI2 (10.31) 



We now turn to ||9^||4. 
\\dS^ || 4 

< c ( gjgi + QoQl + g|g? ii^iIhi 

+ C ( ||fa 2 + eg)0ife||4 + ll^o^H 4 + Q|||^|| 2 (H770II00 + Halloo) 

+ C(g|||7 7 o||oo||97 76 || 2 + 11(770 + e )77 6 %||4 + ||(^ + 9eoK 2 ||4 + || % 2 ^||4 

(10.32) 

To further estimate ||<9^( fe )||4 we use the following estimates of individual terms: 

||ijg0ifr||| < Halloo IMII H\\m 
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|| e§o»?7fo || 4 < ||?7o||~ ||eo|| 2 2 IMItf 1 

i i 

WvodVoVbW* < \\dr)o\\oc IMI«> WmM \\Vbh 

i i 

||Wfc^fe||| < IMU |M|oo WvoU \\Vb\\m 

i i 

||e ^^ 6 ||4 < Halloo ||e ||oo ||eo||| IMItf 1 

i 

H^o^blU < \\r]b\\oo ||&7o||& H^olU \\Vbh 



' 3 
\ 

|<9e ^||| < ||e ||ifi |M|£> WvbU 



3 1. 

2 



\\Vb d Vb\\i < ||%||oo ||%||2 1Mb (10.33) 



Recall that ||?7&(i) ||hi can be estimated in terms of ||eo(i)||.ffi; see (10.8). 

Since rjb is driven by the bound state amplitudes (Qo and Qi), which have dif- 
ferent behavior on the intervals Ij, we now estimate r)b(t) separately on I = [0, to] , h = [to, t{\ 
and I2 = [h, 00). 

We now introduce appropriate norms on different time scales. Define 

3 

M Io (t) = sup ||e (0l|oo + sup (t) \\r/b(t)\\oo + sup ||e (t)||#i 

0<tAt 0<tAt 0<tAt 

+ sup ||9e (*)|| L? + sup (t) 2 \Q (t)\ + sup iQ^t)] (10.34) 

0<tAt 0<tAt 0<tAt 

3 

M h {t) = sup (t)2 ||e (t)||oo + sup IMOH00 + sup ||e (t)||^i 

*0<tAti * <tAti t <*Ati 

+ sup (t)i ||de (*)IL« + sup |Qo(*)l + SU P ( 10 - 35 ) 

*0<tAti t <*Ati t <tAti 

M /2 (t) = sup (t)§ ||e (0l|oo + sup (*-*i> 2 ||?7 6 (OI|oo + sup ||e (t)||i/i 

*1<*AT ti<tAT ti<tAT 

+ sup (t)§ ||9eo(0llL? + sup |Q (t)| 

+ sup (i-ix) r'|Qo(*i)| |Qi(*i)| (10-36) 

ti<tAT 



Remark 10.1. By Propositions 10.3, 10.4 and 10.5, the e contributions to the 
norms Mj k are controlled in terms of the initial conditions. Therefore, to control 
M Ik , it suffices to hound Qj and r/ b . 

The above estimates can be used together with the bounds on Qj of section 9 to 
obtain the following three propositions, which give bounds for r/b on the intervals 
J , 1 1 and h- 
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Proposition 10.7. (rj b (t) for t G Io) Assume t G Iq — [0,t ], i-e- <5o(0 < 
£ (Mx + £ ) (t)~ 2 for i G [0,* ]. Then, for ||0(O)|| X sufficiently small 

IMOIIoc < C( ||0(O)||x,M /o (t o ) ) (t)- 1 (10.37) 

Proposition 10.8. fr? 6 (t) fort G /J Assume t e h = [t ,ti]. Then, for \\(f>(0)\\x 
sufficiently small 

H(t)\\oo < C( \\m\\x, MrM^^h) ) (10.38) 

Proposition 10.9. (rj b (t) for t G h) Assume t G h = [h,T\. Then, for \\</>(0)\\x 
sufficiently small 

IMf)IU < CilimhiM^M^M^T)) (t-fi)"* (10.39) 

In our estimates of section 11, we shall use the following result to estimate 
the size of correction terms in the system for Q and Qi for t > t , where Q is 
monotonically increasing. 

Proposition 10.10. Let ( = ((x,t), x G 1R 3 , t > t , with ((x,t ) = satisfy 
the following dispersive equation 

idtC = Ho( + P c ( s 1 (t) + s 2 (t)( + xC 2 + C 3 ) , (10.40) 

where for all k > and j = 1,2: 

= o(c([0(o)]xQ o 5 (O 

<9 t Qo > for t > t , and Q < £ . (10.41) 

Suppose HCWHffi < WWl^M^^M^T) ) for all t < t < T, where 
[0(O)] X is sufficiently small. Then, ||C(*)||oo < C ( [0(O)] X , M 7l(tl) , M h (T) ) Q (t) 1/2 . 

Corollary 10.1. Let p > 6. £™ = sup t , <t ei(t'), then 



Corollary 10.2. Let p > 6. 

which follows by using 

Wxe- lHt PcX9h < c(t)- 3/2 \\9h. 
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Corollary 10.3. 

\\d\l\oo <CQl/\t)^^^ 2 ds. 

Proof: This result follows from applying d k to the equation for ( and estimating, 
as above, in L p for any a and p > 6. By the Sobolev inequality this implies control 
derivatives in L°°. 



Remark 10.2. In our applications of Proposition 10.10 and its corollaries, ei(s) 
will be given by the source terms depending on Q and Q'(. For t > t±, Q\ = 
0({t — Therefore, since the lowest order term in Qi contributing to e 1 (s) 

is 0(Qi) ( see equation (7.32) ), it follows that for t > ti 

H\\w^ = 0((t - + 0((t - h)- 1 ) + 0((t - tx>-3). (10.42) 

Since rj = O(r] b ) + O(r] ), the conclusion of the main theorem for large t, t > t 1} 
follows. Namely, for 

|M| w *,oo = 0((t - ii>-3) + e- inot P c [r](Q) + tfl (10.43) 
where the non-free wave part is coming from spatially localized source terms. 

11. Beginning of proof of Proposition 9.1 

The key to Proposition 9.1 is the following more detailed version of Corollary 7.1: 

Proposition 11.1. Let t > 0. The equations for P ,P\ can be written in the 
following form 

^ = 2TP P^ + R^[r ]0 ] + R^[r ]b ] + R^[P ,P 1 ] 

^ = ^rPo^ + ^M+^M+^^o^i]- (11.1) 

where 

(i) Rq^ [rjo] are 7] -dependent terms only both local and nonlocal in time, and may 
also depend on ao,Pi- 

(ii) rj b is the bound state driven part of the dispersion; see section 5. 
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(in) R 2 depends only on Pq, Pi, t, but not on r/; It is formally linear in Pq, of high 
order in P 1; and contains both local and nonlocal in time terms, 
(iv) depends only on P , P\,t, but not on rj. It is of high order in P 1} local 
and nonlocal terms included, but has terms which are linear in Qf ~ y/Po- 

The proof uses repeated application of near-identity transformations of the 
variables <So,/3i, derivable by integration by parts ( see the discussion of resonant 
and nonresonant terms in section 7 ) and the decomposition of r\: i] = i]q + eo + r/b; 
see equation (10.1). 

The proof is long so we break it up into three parts, which are presented in 
three different sections. The following is an overview. 

Part 1: The terms i? °^[?7o] and R^\f]o] are forcing terms in the ODE dynamics, 
which are driven by the dispersive part of the initial conditions. They are studied 
and estimated in this section; Proposition 11.2 

Part 2: The terms R^fa] and i?^^] are studied and estimated in section 13; 
Proposition 13.1. 

Part 3: The terms ^ 0) [P , Pi] and ^(Pqj-Pi) 

are studied and estimated in 

section 12; Proposition 12.1 

In Parts 1, 2 and 3, we require estimates for allt > 0. On I = {t : < t < t } 
we use the d priori bound on Po{t), implied by the definition of to; see eqn. (8.5). 
For t > to, we use the 
Monotonicity Property Q 

Qo and ^ are montonically increasing, (H-2) 
Qi 

where Q and Qi are the modified bound state energies related to P and Pi ( see 
(9.7), (9.1), (9.2) ). This monotonicity property shown to hold at t = t and is then 
shown to continue for all time, t, by a continuity argument; see section 14. Since 
there are many terms, we focus on those which are most problematic, namely, 
those which are nonlocal in time and of slowest time-decay rate. This calculations 
are very lengthy and before embarking on them we present a calculation, related 
to the normal form discussion in section 7, and which is repeatedly in order to 
exploit rapid oscillations in time. 

Expansion of oscillatory integrals, resonances and improved time decay 
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In deriving and estimating the terms in (11.1), we must frequently expand 
and/or estimate terms of the form: 



X 



1: J Q e- mo ^P cX 2(s)e i ^ s ds S j (11.3) 



Here, X i an d X2(s) are localized functions of x, X i is independent of s and X2{s) 
depends on s through its dependence on a (s), /3i(s) or 770(5). 
Recall that H , define in (4.6) is of the form 

Ho = Ho + E a 3 (11.4) 

and Hq = o" 3 (— A) plus a matrix potential which decays to zero rapidly as \x\ 
tends to infinity. In (11.3) we would like to integrate by parts, exploiting the 
oscillation of frequency E . However, "peeling off" these oscillations is a little 
tricky because cr 3 does not commute with Hq. We handle this as follows. 
In [10, 11], it is shown that 

ZH W = a 3 (-A - E ) , 
Ze- iHotPc W = e -W-A-Eo)t ( 1L5 ) 

where W denotes the wave operator 

W = lim e -in t e -ia 3 (-A-E )t 

Z = W' 1 (11.6) 
By (11.5), we can rewrite (11.3) as 

'xi , e- iHot J*We mi{ - A - Eo)s ZP c e lUlS X2 (s) ds^ 
' X i ,e- iHot J*We- ia3 (- A - Eo >ZP c e lUlS X2 (s) ds'j 
' X i ,e- lHot J*We ia3As Z ■ We ia3Eos ZP c e lQlS x 2 (s) ds 
' X\ , e ~ iHot f* We iasAs Z ■ We i{a3Eo+mi)s ZP cX 2{s) ds 

Xi ,e~ iHot fwe^'Z ■ {-i)^W{<j z E Q + m l )- 1 e i{a3Eo+m ' )s ZP cX2 {s) ds 
K Jo ds 

'xi ,e~ lHot f Wa^e la ' iAs Z ■ W(a 3 E + m 1 )- 1 e i ^ Eo+m ^'ZP c x 2 (s) ds 
\ Jo 

'xi , J* e m ^H ((T 3 E + mj-^'Pc X 2(s) ds^ (11.7) 
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In the previous string of equations we have used the notation / ~ g to mean 
equality up to terms which are local in time. Note that g$Eq + IVt\ is invertible 
since its determinant is Vl\ + Eq. We can therefore carry out this procedure any 
finite number of times to arrange, up to local in time terms, an expression involves 
an operator of the form x exp(iH,o(t— s))HqX> k > 1, where x is spatially localized. 
Therefore, the enhanced local decay estimate (4.43) of Theorem 4.3 applies. We 
shall use these observations, together with the detailed dependence of X2{s) on on 
a , (3\ etc. to control certain terms in 



Part 1: Estimation of 



4% 



and 



Proposition 11.2. Assume either t < t or Montonicity Property Q on 

[t ,t + 5*]. Then, 

4 0) M|, |4 1} M| < ^(t o ,Nx)(t)- 2 + 0(^o)2rP o P 1 2 (11.8) 
bo(to,[Vo\x) = O((t ) _1 Mx) 
M*o,fa>U) = Wo>^Mf) 
Proof of Proposition 11.2: 

Proof of esimate (11.8) for [r] ] : The key terms are those in the Po equation, 
which are decaying most slowly with t. These are linear in r] , since ||^o||oo = 
(9(t~ 3 / 2 ). We focus on the most difficult terms. These are nonlocal in r/ (t). Recall 
equation (7.52) for a and that the equation for P , (7.6), is derived from the <So 
equation (related to (7.52) by a near-identity change of variables) by multiplication 
by a and taking the imaginary part of the d t a equation. 

We consider the following representative "most problematic" terms in R§\r)o], 
whose estimation introduces the necessary methods for treating them all: 

(Tl) O(a ) ^,/ o < e-^''C- s )p c M/ ( S )^ 1 ( S )f7 ( s )rf^ O (o^e^ + |&| 2 ) , 

(11.9) 



( also with ^f ^in replaced by ^o^i^o, ^o^i^o ) 
(T2) 0(ao)0(aoPi) (x, [ e^'^P^o) ds. 



(11.10) 



Estimation of Tl: Since the time-integral is bounded by O(£o)(t) ' [Vo]x, by 
the Cauchy-Schwarz inequality the second term in (11.9) is bounded as follows: 



0{a,)0{\h\ 2 ) (x, jT 



<o{£ Q )p Q Pt + o{£ ){t)-*\\r, \\ 2 x . 



;ii.ii) 
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We now control the first term in (11.9). O(aoj3i)0(ao). We argue that the 
key contribution from this term which must be bounded is of the form: 

OiaoWoPie"^^ (x, f Q e~ m ^ P^^^^ds^j (11.12) 

To see this, consider the term in the a equation which corresponds to the first 
term in (11.9): 

(x,^e- Wo(t - s) P c *o(s)*i(^ ( S )^)c»(^ie- iA +*), (11.13) 

Next we integrate with respect to t, and integrate by parts, making using the 
oscillatory exponential factor. The result is a boundary term, which can be sub- 
sumed in the definition of ao, by a near identity transformation, followed by a 
time-integral to t. The latter contributes terms to the Pq equation (which has 
been modified due to the slight redefinition of ao) of the following type: 

O(a o )O0 1 d t a o + «„$&) (x, f e'^^P^^ds^ (11.14) 

O(a )O(& Pi)^- t (x,/ o *e-^ o( *- s) P c * *i%^) (11.15) 

Since d t a = C(|a | 2 |/9i|) and dtfc = O(\a \ |/?i| 2 ) (11.14) is bounded by 

0(80)0(^(1}-! [ m ] x ), (H.16) 

where we have used 

Ikolloo < CMx (t)' V2 - (11-17) 

By the Cauchy-Schwarz inequality this is bounded by 0(So)(2TP P^ + [vo\x(~t)~ 3 )- 
Therefore, the contribution from this first term satisfies the estimate (11.8). 

Obtaining a bound on (11.15) is more involved. We use the local decay estimate 
of Theorem 4.3: 

\\ X e- mot P c n f\\ 2 < (t)-" 2 \\d k (xYfh. (11.18) 

The expression (11.15) bounded by: 

= 0(a Q )0(WM (<X, *o*i%> - i(x, [ e- mo ^P c H y yirjods)) 
< O(£ ) ( TQoQj + (ty 3 ) 

ft rla 1 

+ O(a 2 A) [77o]x jf (t _ s)5/2 ^||5o||AI- (11-19) 
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The latter integral requires detailed estimation on different time scales. 

To estimate the last integral, we split the range of integration into three re- 
gions: 

h = {s : < s < t } 
h = {s:t <s< h} 

J 2 = {s : ti < s < T} (11.20) 

Estimate on J : Assume t < t . Recall that by (8.5), |«o(s)| < 0{Sq){s)~ 1 . 
Using this we have 

~ f* 1 1 

0{6i 2 M lno]x J o ( t _^ 5 /2 \ a °\ ds 

= Oda.1*) (4ra|A) (M? /J ¥ _L_-i T | ao|/ 3 1 | ds ) 
< o(sl) [2rp„p, 2 + W\f(t)-i-i\ , 

where we have used that ah < a 4 /4 + 3fet/4. 

Estimate on 1^. Let t be such that to < t < t\. We break the integral into an 
integral over [0,to] plus an integral over [to,t\. Recall the definitions of Qj(t), j = 
0, 1 in terms of \&j(t) displayed in (9.7). Using that 

• for s G I , \a {s)\ < O{S )(sy l and for 

• s e h, Q is increasing and Qo(s) < £oQi(s), 
we have 

0(\a (t)\ 2 \m\)[Vo]x ( /%£ ) (t _ s)5 1 /2(s) 3 /2 IA(^)l ds 

< I type bound + O(\a \ 2 |&|) [r)o]x £ ^r^j^i ds 

< J type bound + C(|a | 3 |A| PilU) Mx jf ^ _ ^5/2^3/2 ds 



= I type bound + O (^£ 2 KI • Kl • KHAIJ [i]o]x (t) 2 ds 
= J type bound + O (s ■ Ql ■ £ 5 Qi) [r)o]x (t)^ 
= lo type bound + O (s^) Q0Q1 



= J type bound + O ( QoQl + (t)~ 3 ) , (11-21) 
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which is a bound of the type in (11.8). 
Estimate on I 2 : 



to 



"oil (k 

s )3/l 

fiollA 



= I & h type bounds + <9(agA(f))fook £ ^ _ ^5/2^)3/2 rfs 

(11.22) 

The latter integral must be treated differently from the previous terms. For 
t e I\i we used that Q monotonically increasing and bounded by a small constant 
times Q\ to treat terms perturbatively. On J 2 , Qo dominates Q\ (which decays) 
and we must use a different argument. We return to expression from which the 
last term in (11.22) is derived: 

0(a o )0(3^i) ((x^o*!*)-^^^^""'^^*^*)) • (H-23) 
We need to expand and estimate time integral: 

[\- iHo{t - s) P c HoMs)Pi(s)e- iX+s r] (s)ds = e~ mot f ' e l{n ^ s 'P c H ao(s)Pi(s)rj (s)ds 
Jti Jti 

(11.24) 

which we do using integration by parts. We carry this out, then take the inner 
product of the result with a localized function, x, and then finally multiply by 
C(|a | 2 IAD- The result is of order 



«o| 3 \Pi\ 2 IM*)IU + l«o| 3 (x,e- mot HoX 



-iHot' 



O(\& \ 2 IAD j ti e- ino{t - s) Pc{n,-\ + r l n oX e- lX+s j- s {&o{s)Hs)%{s^ ds 

(11.25) 



The first two terms in (11.25) are bounded as follows: 

|«o| 3 IMOII00 =0(£o) (|«o| |A| 2 ) ll^oWllc 

< O(S ) (QoQj + \\m(t)\\lo) 

< O(£o) ( QoQi + (t)- 3 ) 

72 



Wlftl (x,e- iHot Hox) < C|« | 5/2 (|« | 1/2 IAD (t)- 5 / 2 

< 0{Sl /2 ) (QoQl + (f)-H ) < 0(^/ 2 ) (Q,Q\ + (*>-3 ) . (n.26) 

We now turn to the nonlocal term, in (11.25), which we denote T\. 

Jx ~ O(\d \ 2 IAD fe'^it-s) Pc{Ho _ X+ )- l Ho X e- tX+s 

Jti 

( 9 s ao(s)A(s)?7 (s) + <5 (s)<9 s /3i(s)?7o(s) + a (s)/3i(s)<9 s ?7 (s) ) . 

(11.27) 



Using that 



d t a ~ 0(#a o ), $A ~ O(doA), and <9 t r/ 



-iH r] , 



(11.28) 



we have 



J! ~ O(|a r IAD 



Jti 

( «oA 3 r/ + aoA 2 %( s ) + M s )A(s)7io??o( s ) ) 
= J la + X lb + Jic (11.29) 

Each of the three terms Xy, j — a,b,c satisfies a bound of the form: 

Pal < CEl [770(0)]^ y^yy^p ( 1L3 °) 

We illustrate this by estimating X\ a \ the other two terms are estimated simi- 
larly. 

1 

Using that I a 1) IAI — O(Sq) we have 

rt 1 



|Zla| < l«0| 2 IAI 



ti (t- s)a 



|A| 3 |a | btol 



< £ 2 MO)]^ (sup(*-^|A(*)|) /* . 1 ... 1 .3 7^ 

\t>tl J Jti (t - S)2 [S - tl)* {S) 



Separate estimation of the contributions from the intervals [ti, and [\{t- 

ti),t] yields the bound: 



\Iia\ < C El [ m (0)] x sup(*-^|A(i)| 



t>tl 



1 1 



(11.31) 
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Estimation of T2: Consider the term 

T2 = O(a o (t))0(«o(t)A(t))(x, j ( < e^('- s )p c% 2 ( S )r ?o ( S )^ ds. 



(11.32) 



te h = [tojh]: For t G h, 



|«o| 2 |A| < c£ 2 |a | |/3i| 2 ; Proposition 9.4. 



(11.33) 



Therefore, 



|T2| < c£ 5 |fi | |A 



Jo 



< c£* lorol |A| 2 <t) _3/2 sup( S )i||?7 2 ( S )77o(s)|| W MnLi 

\0<s<t , 

< cf| |a | IAI 2 (*> _3/2 IMIk 2 Mx 

< cSSimU |«b| |&| 2 (t)- 3/2 

< <[^ ]x ( PoA 2 + (t>- 3 ) 



t E I 2 = [h, oo): 

To see the relevant terms for £ > t 1; we integrate by parts and obtain, besides 
easily estimable local terms and terms with faster time decay, 

0(«oA) ( X, f e- m ^- s) H PcV 2 b Vo ds^ t>h. (11.34) 
Consider the contribution to the integral in (11.34) coming from s G [0,t ]: 

O (ag(t)A(f)) (x, e- ino ^H P c rj 2 b rj ds^ t > h > t (11.35) 
Consider, the inner produce 

R(t) = (x,j* e- in °^H PcV 2 b Vo ds^ (11.36) 

We have 

|(11.35)| < C |a„(t)| 2 |A(*)l \R(t)\ ~ CQ (f)Q?(f)W)l (U-37) 
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and therefore it suffices to prove that 

\R(t)\ < CQUt) (11.38) 
To prove (11.38), recall that for t > ti y 

' TQoQl 



dQo ^ 2 



eft 

2rQ Q?, (11-39) 



dQi ^ 2 



and 



Therefore, 



dt 



£oQi(h) = Qo(h). (n.40) 



Q M x i I ?fff w (1L41) 

1 + r/ ti Q (s)ds 



Therefore, to establish (11.38) we need: 



2Qi(ti) 

i + rjlQ ( s )ds 



W)\ < C [ - ; I (11-42) 



We consider two cases 

Case 1: T sup se [ tljt ] |<5o(s)| \t — h\ < 1, and 

Case 2: T sup se [ tl t ] |(?o(s)| \t — h\ > 1, where it suffices to prove 



m)\ < (JnrAr ) ("•«) 

In case 2 we prove the bound on R(t), (11.38), while in case 1 we prove that the 

3 3 

expression (11.34) of order (t — ti)~2 (t — 1 )~2 . We first handle case 2, the bound 
(11.43). From (11.39) we have 

j t ( Q + 2Q 1 ) x 0, t > h 
Q (t) + 2Q 1 (t) x Q (*i) + 2Qi(ii) 

= (2 + £ r )Qi(ti) (11.44) 



Therefore, since 



1 + FJ ti Qo{s)ds 
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we have, as t — > oo 

Q (t) x (2 + ^)g 1 (t 1 ) - 2Q.it) 

— (2 + 5 r )Qi(ti) (11.46) 

Therefore, for £ > ii {t — t\ large enough) 

(l + £ r )Qi(*i) < Q (*) < (2 + ^)gi(t!) (11.47) 

and therefore 

Qi(ti) > Qi(ti) 



1 + rf t \Qo(s)ds - i + r(2 + £j) |t — ti| Qi(*0 

- 3r |* — t x | _ 3r \t-t \ v ; 



W)\ < C(x,[°- r3\\HoV 2 b (s)vo(s)\\ L i ds 

\ JO t-s)i 



On the other hand, 

r*o 1 
(t-s 

fto 1 

< C sup H(t)\\ H 2 [r)o\x / rr—jds 

re[0,t ] ^0 (t_ s ) 2 ( S )2 

< CgoMx . 1 ., (11.49) 

{t-t )i 

The bounds (11.49) and (11.48) imply (11.43). 

We now turn to case 1. In this case, rsup s6 [ tl t ] Qo(s) \t — t±\ < 1, 

Q.W < 1^ (n.50) 

and therefore by (11.34) and (11.50) 

|(11.34)| <£ 2 Mx sup ((s-ti)*QiO0) 1 1 1 . (11-51) 

We begin by noting that by Proposition 10.10 and its corollaries, for p > 6 

w*.* < O(£ )(t)-\ t<t 

< 0(£S)Qo(t) 1/2 , t <t< tl 
= 0({t-t!)-*),t > h (11.52) 
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For t < ti, the previous arguments with the known estimates on rjb, and the 
facts that |<5o(0l — ^(^o)/ (t) on I and Qo(t) < SqQi on l\ imply the necessary 
bounds. Collecting all these, we have 

rt 



( x , jW-^PcHovlmds) <O(S ) 2rQ Q 2 1 + (t )- 1 (ty 



Proof of estimate (11.8) for R^: The key terms to consider in the Pi equation 
are the slowest decaying nonlocal terms; see (7.39) and (7.32). Since b\ enters as 
b\ in (9.22), we need to bound R { 1] by O(S ) [{t }- 1/2 {ty 2 + 2rP A 2 ] for t > t . 
The slowest decaying nonlocal term in the Pi equation arises from the balance: 

idtPi ~ Mipu,^ 2 )\Pi\ 2 e iX+{T)t (11.53) 

Since the Pi equation is obtain by multiplying by Pi and taking the imaginary 
part, we must estimate: 



11.54) 



the leading order nonlocal part of which is 



Pi\ 2 Pi(x, / e- ma(t ~ s) P c ^ ^iri ds)e lX+ 



\Pi\ 2 Pi(x, / e-^'('-)p c «oAe- lA + s ^o(s)^)e iA +* 



Integrating by parts (using the oscillatory factor e lX+t ) and removing nonresonant 
local in terms by near identity transformations the key term is 

ft 



\Pi\ 2 Pi (x,f Q e^^HoPcaohxe^^dsy^. 



(11.55) 



Suppose t > t\. Since we have no factor of q.q outside the integral (local in 
time ao(t) factor), the estimate for I = {0 < s < t } is the critical part and 
requires the 0({t)~ r ) bound on I and Theorem 4.3: 

^ m ^n P c a PixVoe- iX+s ds X 



X- 



< 



to 



ds 



(t- S ) 5 /2( s )3/2 

*o ds 



\a (s)\\Pi(s)\ + 

+ [Vo]x f 

J t 



/" + /' 

J to J t\ 



• • • ds 



C[vo]x f 
Jo 

< O(£ )[r) ] x Jo {t _ s)5/2{s){s) 3/2 < lvu JA JtQ {t _ s)5/2{s)3/2 

< O(£ ) f(t)- 5/2 Nx + 0(^o)Nx|5o(t)|(t)- 3/2 + Mx|5 (t)| (*-*i>-'<*>-5 



ds 



\a (h)\ IA00I + 
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Note that we can extract the factor \ct(t)\ from the integral for s > to since in 
this range Qo(t) is monotonically increasing and Qo(t) ~ Po(t) ~ l^oWl with 
correction terms which are which are rapidly decaying in time and which are 
therefore dominated by the first term. Multiplication by (9(|/3i| 3 ) prefactor gives 
the bound 



C(£ , IVoU) 



11.56) 



This completes the proof of Proposition 11.2. 



12. Local and nonlocal ODE terms: {Pq,P\) of Proposi- 
tion 11.1 

In this section we prove estimates on the terms R^\Pq, -Pi) an d R^\Po, P\) of 
Proposition 11.1. 

Proposition 12.1. Assume either t < to or Monotonicity Property Q on [to, to + 
5*] . Then, For m > 4, 

(1) |4 0) (Po,Pi)| < O(SS)TP P* + 0((t)- 3 ) (P ,P 1 « 1) 

(2) |4Vo,Pi)l < O(£S)TP P^ + O(a )P^ + O((t)- 3 ), \a \ ~ 

12.1. Proof of part (1) of Proposition 12.1: 

The most problematic terms are nonlocal, slowest decaying. The terms which are 
linear in r\ in the a equation contribute the slowest terms; these are nonlocal in 
time, t. 

We have to consider terms arising in equation (7.52) of the type: 

( X ,v)e- tX+t a-o(3 1 , (x,v)\Pi\ 2 , (x, v) fte' 2 ^ , (12.1) 

where 2 ~ V and n(t) = r) (t) + e (t) + i] b (t); see (10.1). 

As calculated earlier, the last term is resonant and its contribution is +2rP -Pi 2 
to the P equation; see (7.55). 

The leading ODE terms in n are the source terms of the type in (7.32): 

|* | 2 *i,*o^*o|*i| 2 ,^?; (12.2) 
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Recall that the Pq equation is obtained from the «o equation by multiplication by 
aio and taking the imaginary part. 

Solving for 775 and plugging the source term contributions into the (%, 77) terms 
in the Pq equation gives, apart from the resonant term, terms of the type 

mO(a Pie iX+t + \f3i\ 2 ){x, f e- iHo{t - s) P c ( |^ | 2 ^i + *o*i + ^o|^i| 2 + ? ) ds 

(12.3) 

For t > ti the slowest terms are C ( | \I/o 1 2 ^I/i + ^o^i) source terms. For to < t < ti 
the problematic terms are C(\I/o I 1 2 + ^0^1) source terms, since on this interval 
f3\ does not necessarily decay. Since ^ ~ e~ lX+t f3\ + higher order terms, we can 
integrate by parts the 0(|\I/o| 2, I'i + ^0^1) terms in (12.3) to get arbitrary order 
/3™ terms, which are nonlocal. The remaining terms are local and its lowest order 
term is of the order 

a O(a e^ t p 1 + |A| 2 ) ( X , 0(K| 2 /3ie a+ ' + a^e'^)) . (12.4) 

The nonresonant local terms can be transformed to higher order by a near identity 
change of variables giving 

O(\a \ 3 Pi m e int ), tt > and m arbitrary large, 

while resonant terms are of the type already derived but of higher order (bounded 
by 0(£o)QoQl)). We are therefore left to consider the following nonlocal terms 

a O (a /3ie- iA+t + |A| 2 ) (x, f Q e~ m ^ P c ao(3\ m ^ s x^s) (12.5) 

Consider the first term in (12.5). Due to the oscillatory factor e~ lX+t in the 
0{aol3\e~ lX+t ) term we can, by a near-identity transformation, in the ao equation, 
remove this term in exchange for one higher order 

+ higher order + local terms, 

and by further near identity transformations, we are left with the following slowest 
term 

0(K| 2 /Wr (x, Jl e- m ^ P c ... ds^. (12.6) 
All other terms are of order o(£ )O(P Pf), o(£ ){t)~ 3 or higher order, for £ small. 
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dt(x, f e- mo{t - s) P c a^l m xds) 
Jo 

X ,a (t)f3l m (t)x) ~ (x, i[ e- in °^HoP c xaof3l m e ins ds 



The first term is local and will contribute o(1)0(PqP 2 ). It remains to estimate 
the nonlocal term. Note that 

/' e- m "^e iujs xa pf m ds = e~ m t e~ t{HQ ^ ){t ~ s) e l{u)+1 * )s xa f3 2m ds 
Jo Jo 

Consider now the second term in (12.5). Using the oscillatory factor e~ luJt , one 
can transform the O(a \f3i\ 2 )(x, I • ' •) to higher order. 

We now turn to the second term in (12.5). First, let's consider the case where 
t > t x . 

O(a \Pi\ 2 )d? ( X,J*e- in °^xa f3l m ds } (12.7) 
We consider (12.7) for t>t\. For this we require the following 
Lemma 12.1. For t > t±, 

Qi(t) > Qi(h) [1 + 4(r' + 5)Q 1 (t 1 )Q (t 1 )(t - h)}- 1 (12.8) 



Proof of Lemma 12.1: For t > ti 

dQ 
dt 

Therefore, 



l - > -4(r' + 5)Q Qi (12.9) 



dQi 1 > -A{Y'+5)Q {t) > -A(V+6)Q {t 1 ) or ^ < 4(r'+5)Q (t 1 ) (12.10) 



dt ' - ,•// 

since Qo(ti) < Qo(t) (Qq | for t > t ). Integrating from t to t±, we get 

Qiitr 1 - Qiih)- 1 < 4(r' + 5)Q (t 1 )(t - to, (i2.il) 

which is equivalent (12.8). This completes the proof of Lemma 12.1. 
Estimation of (12.7) for t > t\ Carrying out the differentiation in (12.7) we 
find that it suffices to bound terms of the type: 

O(a \p x \ 2 ) (x,P cX ) «oA 2m = O(£ )Q QlO (a |A| 2 ) (x, jT e- iH ^P c H^x^l m ds 

(12.12) 
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We now consider (12.12), which we break into the sum of three integrals: 
O(a \(3i\ 2 ) ( X, ( fj + jT 1 + £) e^**-* P c ^oft 2 ™ ) (12-13) 



J* t * : By the local decay estimate for e lHot HoP c of Theorem 4.3, the integral in 
(12.13) is bounded above by: 

{(dV^? 1 (12 - 14) 

This in turn is bounded above by 

Qo(t) 1/2 H^^QTis), (12.15) 

since Qo is increasing for t > to- 

We now aim to further bound (12.15) by "extracting" powers of Q±(t) from 
under the integral. Recall that for s > ti 

Qi{s) < 



l + 2TQ (t 1 )Q 1 (t 1 )\s-t 1 \ 

Q 1 (t 1 )Q (t 1 )2F 



2TQ (t 1 ) 



l + 2TQ (t 1 )Q 1 (t 1 )\s-t 1 



and 



We write 



< miii {Qx^O, '(tiXs-*!)" 1 } (12-16) 

Q- 1 ^) = Qi^i)- 1 ^, and Q x < £„■ (12.17) 



QT(s) = Qr%)Q%) < Qr% h * r , 9rn * r 

(s - ti) fe (2rQ (^i)) fe 

_ i c—rk 

= Q?-~ k (s) 



(s (2rQi(*0)* 
< QT M "(*)V*(s-ti) 4 - (12-18) 

where we have used that Qi(^i) > Qi( s )- 

We consider separately the cases t > 2t\ and ti < t < 2t\. 
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t > 2ti m . Take k = 2 and m > 2(r + 2). Then, the integral in (12.15) is bounded 
by 

£o Qo(t)Ht)~ min ^\ p>0. (12.19) 
This implies, for k sufficiently large, 

0(«o|A| 2 )(x, jf •••>!< C(^o) ( QoQ? + ) • P > (12.20) 



ti <t < 2t 1 : Let t = ti + 2M, M = {t - ti)/2 and rewrite the integral as 

/■* / fti+M r t!+2M \ ,7 S 

In a manner similar to the previous estimate, using that Qi is decreasing we have, 
by (12.18): 

rti+M rl q , _ 

/ 7T^uQi(sr /2 (s-ti)- fc < c(f-(f 1 + M))-*Q 1 (f 1 ) m/2 

< ^-^-^(^gi^f- 1 

< c(t - t!)- fe [ 1 + 4(r' + 5)Q 1 (t 1 )Q Q (t 1 )\t Qi(*)Qi(*i) m/2-1 

< 0{£ Q )Q 1 {t) 

for > 1. Hence 

Oia^^Ql' 2 / • • • < 0(£ o )QoQ 2 . (12.22) 
Furthermore, by (12.18) the integral over \t\ + M : t 1 + 2M\ is bounded by 

r t 1 +2M , , , 

/ „ - — u ; - ^ i(*r /2 < c gi(*i)-/ 2 (M)- fe < cQ^r^t - h)- k 

Jt!+M (t — S) k {S — ti) k 

(12.23) 

and, as above, using the upper bound (12.8) for Qi{t\) in terms of Qi(t) we have 

rti+2M 

/ ■ ■ ■ < O(S )Qi(t) (12.24) 

for fc > 1. Therefore, for all t > t\ the nonlocal (and local) ODE terms in the Qo 
equation are bounded by 

0(£o) [QoQ 2 + (t>~ 3 ] , (12.25) 
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provided we control the integral /q 1 • • • ds. 
Jq 1 : Consider 

lit) = Q 1 /2 Q 1 £ 1 -^- k Q 1 / 2 Q^Q 1 ds. (12.26) 

By the mean value theorem, Qi(s) = Qi(s) — Qi(t) + Qi(t) = Qi(s)\t — s\ +Qi(t), 
where s < s < t. Then, 

l rti Aa l 

lit) < QUt)Ql(t)J o -^-^- k QUs)Q?-\s) ds 

+ h(t), where (12.27) 

1 rh (7s I I 

h(t) = cQUt)Qi(t)J Q _ s)k _ 1 Ql(s)Qr 1 (s) QS(s)Qi(s), (12.28) 
where we have used 

Q 1 = o(Q%Q 1 > )+h.o.t.. (12.29) 
i 

If s < t , then using that Q§(s) < S^s)^ 1 < Sq(s)" 1 , we have the bound 

h(t) < O(SS)Ql(t)Q 1 (t)(t}- 2 <O(SS)(Q (t)Q 2 1 (t) + (t}- 4 ) (12.30) 
If s > t , then since Qo(s) is monotonically increasing for s > t , and we have 

hit) < /„" J^QlWQT-Hs) Q lW (12.31) 

We now expand the latter factor of Qi in the integrand using the mean value 
theorem. Specifically, there exists s' with t < s < s' < ti such that: 

Qi(s) = Qi(fi) + gi(s')(s'-fi) 

= Qi(*i) + o(gl(s')Qi(s'))|s'-*il 

< Qi(ti) + o(Ql(*i)Qi(s'))|s-*il, ( 12 - 32 ) 

where the last inequality follows by monotonicity of Qq. 
Substitution into integrand in (12.32) gives: 

h(t) < cSS ^f^- QUti) jT 1 ¥ J ^z- 3 Qlis)Qr 1 is), (12.33) 

where we have used (12.17) to replace Qo(ti) by Qi(ti). A higher order term (one 
proportional to Qi(ti)) is subsumed by the constant, c. 
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From (12.8) we have 

Qi(h) < (l + A(T f + 5)Q 1 (t 1 )Q (t 1 )(t-t 1 )) Q l {t) 



12.34) 



Case 1: 4(r' + 5)Qi(ii)Q (*i)(*-*i) < 100 Here, Qi(h) < 101Q(f), and therefore 

ds 



I l{ t)<cES Q ° mm " 



ms)Qr\s) 



12.35) 



(t-h) Jo (t-s) k - 3 

We now split the integral in (12.35) as J h = j* *° + j£. Using the (s)- 1 decay of 
i 

Qq for s < t we have 

f° ■ < £§{8)-\ (12.36) 
Jo 

and using the monotonicity of Q for t 1 > s > t and the relation Qo(l) — 
SoQiih) < 101£ r Qi(t) we have 



I* 1 ■ ds < sgQZih) = v / Toi£ p r *Ql{t) 

J to 



12.37) 



Therefore, choosing m sufficiently large, we get 



hit) < c> (s^Qom! mt)- 1 +s p o 2 Qoit)Q\it)) < c"£S 3 (QoQi + (tr 4 ) , 



(12.38) 

where = min{p 1 ,p 2 }- 

Case 2: 4(T' + 5)Q 1 {t 1 )Q {t 1 )(t -t ± ) > 100 Then, (12.8) and monotonicity of Q 
for t >t implies 



(f-fi> 



<A(T' + 5)Q (t)Q 1 (t). 



(12.39) 



This gives 



h(t) < C£SQ (t)Ql(t) (12.40) 
We now turn to (12.7) for t G [<o,*i]- It suffices to estimate the nonlocal terms 

0(«o|/9i| 2 ) (x, jT* e- iW °( t -)p c 7Cx«o/3i 2m ^) (12.41) 

for to < t < t±. In this region Q and (Qo/Qi) are increasing functions. Also 
Qo(0 < £o"Qi(t). The main difficulty is the need to pull a factor of Qi(t) out of 
the nonlocal term. To this end we use the following proposition: 
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Proposition 12.2. There exists a constant 5 > such that for t > t 
2r [ t Q Q?ds<5Ql(t)Q l (t) m - 2 + Qo{t)Q r r\t) + h.o.t. 

Jt 

Corollary 12.1. 

£ j^^qT{s)q?-\s) < c(i + s)Qfar k . 

The corollary follows from Proposition 12.2 by the Holder's inequality. 
Proof of Proposition 12.2: Recall that 

dQ{] 2TQ Qi + R , ^1 < -ATQoQj + O (qIq? 



(12.42) 



(12.43) 



dt 



dt 



Claim: 



-4 [\ds < O(£ °)\l + f Q (s)Ql m (s)ds 

Q2 Jto I Jto 



12.44) 



12.45) 



Proof of claim: The leading order term of Rq, in the variable oto, which is nonlocal, 
is a term of the form (12.5). From this term we have after integration by parts to 
obtain e'^^'^Ho from (12.5), 

fRo < QUt) /'qx(0 f d \ 5/2 Ql /2 Q?dsdt> 

Jto Jto Jto \t S) i 

rt' ds 



< 



l rt rf W s j /-i' 



< cQUt) I Q^dt' (t-t')- M ' + (t-t a )- M ' + (t-t'}- 1 - s ' [ Q a Q\ 

< cQl(t) 
thus proving the claim. 



2 m 



ds 



O(£ )(l+ / QoQi m ds) 



We first rewrite the right hand side of (12.42) and use the differential equation 
for Qo and the above claim to integrate by parts: 



2r f QoQTds = f 2TQ Q\Q™- 2 ds 

Jto Jto 



/*or 2 4[Qo- [ s Ro 

j to as \_ j to 



ds 



= Q?- 2 (t') 



Qo - [ S Ro\ ' - f ^-Q m - 2 (s) \q (s)- [ S Ro 

Jto J to Jto CIS L Jto 



ds 
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= Q?- 2 (t)Q (t) - QT' 2 (t )Q (t ) - Q?- 2 (t) f t R (s)d 

Jt 



+(m 



2) f dt'Q?- l (t')Q (t') 

Jt» 



4TQ (t') + O{QUt')Qr 2 (t') 



+(m-2) / dt'Q™-\t' 

Jto 



4rQ (t>) + o[QUt')Q?- 2 (t') 



to 



R (s)ds 



The first term on the right hand side is of the form we want (the right hand side 
of (12.42). The second term is negative so we can drop it. The third term is 
bounded by O(S )QT' 2 (t)Ql (t) by the claim above, plus 0(Jl Q Q™ ds), which 
is of the above form and can be absorbed by the left hand by smallness of £ , 
where we used t < s < t±, Qo < SqQi- This completes the proof of Proposition 
12.2. 



Using the proposition and its corollary it is easy to bound, for t >t , 



O(a \^\ 2 + K| 2 A) d k t (x,j Q e~ m ^ P c xa fi m ds 
by Q(SS) (Q Q 2 + (t)- 3 ) + h.o.t. 



(12.46) 



It remains to consider t <t . In this case we only know that \a \ < k (t) 1 . 
The first term of (12.5) 

a O (a /3ie- a +*) ( X , J* e- 7 ^*-) P c a f3 2m xds) (12.47) 

is bounded above by 0({t)~ 3 ), due to dispersive estimates and the decay of the 
decay of Qo for t < to- 

In the second term of (12.5), 

|2\ ak 



0(«o|/9i| 2 )^(x, / •••) 



O (Qkt)Qi(t)) (xJ*e^^' ) Pci%X&offi m e ia ' ds} , Q + 



(12.48) 



i 

we need to pull Qo(t)Qi(t) out of the nonlocal (integral) term. 
By dispersive estimates, we have the bound 



O (Qht)Qi(t)) {^pl^WITO ds ( 12 - 49 ) 
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To pull the Q\ term we proceed as earlier. By the mean value theorem, there 
exists s, with s < s < t < t such that 

Qi(s) = Qi(t) + Q 1 (s)-Q 1 (t) = Q 1 (t)-Q 1 (s)(t-s) 
= Qi(t)+o(Qks)Qi(s)yt-s) 

= Q 1 (t) + 0(SS(s)- 1 (t-s)), 
i 

where we have used that Qi ~ Q0Q1 f° r t < to- Therefore, the expression in 
(12.49) satisfies a bound of order 

Ql(t)Ql(t) f \s +k \a (s)\Q?- 1 (s) ds 

+Qkt)Qi(t) jT - _ l s)1 _ +k Qks)o (qI(-s)q7'^)) or\s) d S 

= QUt)Ql(t) f . \ 3+k Qhs)Q?- 1 (s) ds 

JO (t - s)2 +k 

+ £S Q§(t)Qi(t)-^ (12.50) 

where the last term, which is bounded by Sq (QoQl + (£)~ 4 ), is obtained using the 
decay of Qo(s) for s < t . 

It remains to estimate the second to last term in (12.50). Estimating the convo- 
lution, using the (s) _1 decay of Qo(s), we obtain the bound O (^Qo(t)Qi(t)^ (t) -1 
which is not bounded by the desired O ( £q{QqQ\ + (t)~ 3 ) )• 

We will obtain the desired bound by turning to an earlier expression, derivable 
from (12.48). The expression we must consider is 

o(Qh)Q 2 M) (x,^- Wo(t - s) P c 7i: fe x«o(s)A 2m " 2 (s)e^rf s ), n^o (12.51) 

We will show that this term is of order £ p Q (t)Ql(t)(t)~ 3 / 2 , which implies the 
desired bound. We proceed as follows. First, by equation (7.3) of Proposition 7.1 
the equation for a may be rewritten as: 

td t a (t) = (c + ir w )|&(:r)| 4 ao(t) + (c + ir u ) (|&(f)| 4 - |A(T)| 4 ) & (t) + F a 

(12.52) 

Introducing af(t) via the equation 

ao(t) = e- u ^ r ^ T ^ 4t d*(t), (12.53) 
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we have 

id t a*(t) = (c + ir„) (|A(t)| 4 - |A(T)| 4 ) a*{t) + e ^+^)\M^ Fa (12 . 5 4) 
The integral in (12.51) can be written, keeping only leading order terms, as 
' X , J* e- in "^P c H k xai{s)(3l m - 2 {s)e^ s ds^ 
X,£e-^(*- s )p c 7Y fe Xi<5o # ( S )^ 2m - 2 ( S )e^d S + . . .) 
X, fe- m ^P c Hlxi (|A(*)| 4 - |A(T)| 4 ) a*(s)(3 2 1 m ~ 2 (s)e^ s ds 



< 0{£l) [ d " (|A(t)| 4 - IA(T)| 4 ) + . . . (12.55) 
Jo (t - s)2 +k (s) v ' 

It suffices to show 
Proposition 12.3. 

|Pi(t)-Pi(s)| < 0<s<t (12.56) 

(s>3 

For the proof we turn to the Pi equation (7.7): 

^ = -ATPoP 2 + ifc 
R, = R 1 [a J 1 ,r j ,t} = 2%(T 1 Fp). (12.57) 

The key term in S^iFp) is the form p^afe™, ^ 0. We have, after two 
integrations by parts, 

r< _ 



|Pi(t)-Pi(s)| < J s PT(s)a*(s')e ms 'ds' + ... 
= /T(^oV)^' (I - / ^ (/W«oV)) ±e^' ds' + dots 

- (w) + ° (W) + O( ' o) T (s " )l|p(t) - p(s//)l II w ds ' 

= o fe) + ° fifr) + ° O T (s ^ |p(t) - p(s/,)i - (i2 - 58) 

Multiplication by (s)5 and taking the supremum over s < t implies Proposition 
12.3 and therewith the proof of part (a) of Proposition 12.1. 
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12.2. Proof of part (2) of Proposition 12.1 

The proof is similar to that of part (a) but simpler, since we can allow for the 

i 

nonlocal terms to be controlled, in addition, by terms of order QqQ™. 
The leading contributions are again nonlocal, linear rjb contributions: 
From (7.39) we read 

idtPi ~ ■■■ + 2X{ X ,^iVb + ^irh}e lX+t a 
+ Ae iA +*2(x,^i^)a 

+Xe iX+t ( X , ^Irfb + ^i?Vb) + h.o.t. (12.59) 

+ X-\t)(A(t) - A{T))X(t)p x . (12.60) 

The first two terms on the righ hand side of (12.60) are easily seen to be of order 
0(ao)Pi m or O(£q)P P±, by integration by parts over the ODE source terms in 
7]},, as before. 

The third term contributes to the Pi equations, after normal form transfor- 
mations and remaining resonant terms: 



3/V A+ ^(x,|A| 2 / e- in °^P c a \Pi\ 2 xd S ) 
Jo 

and higher order/similar terms. 

The leading term, after integration by parts of the integral term (note that 
i? _1 P c is bounded): 



Jo 

To this end we repeat the argument of part (a) to estimate the above integral by 

QUt)QT~\t) + o(r 3 ) + h. .t.(ODE). 

The main new type of term we need to control comes from the last term on the 
term IZp in (6.35-6.36), coming from the difference A(t) — A(T). 
This term contributes to the Pi equation terms like 

3 ( IAI 2 [«o(0 - "SCO ] 5g(f) ) , 9 [e^i 2 (P (t) - P (T))O(a 2 (T))} . (12.61) 
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The term with phase e lwt can be integrated by parts, and gives higher order 
terms; 

The term without a phase requires the estimate of 

ao(t) 2 — a (T) 2 = 2 f a^^^-ds. 

Jt as 

Using that 

^ = caoffe- 2 ^ + h.o.t. 
as 

we can repeatedly integrate by parts to get 

«o(*) ~ a o( T ) = local terms + f t 0(al)Pf m ds + h.o.t. 

where local terms = O(a Pi) and higher, 
clearly, then 

|A| 2 cto(t) [ local terms + h.o.t.] < O(S ) [P P 2 + (t)~ 3 . 
So, we need to estimate 

Watit) £ O(a 2 )Pl m ds by O(£Z)Q Ql + h.o.t. (12.62) 
For t > ti \/3i\ 2 < Qi is monotonic decreasing. Hence 

| [ T O(a 2 )f3 2m ds\ < |A(t)| 2 [ T QoQ^ds. 
Jt Jt 

For t > ti, Qi I is bounded by 

Qi(t) < Qi(*i)(l + 2rQ„(*i)Qi(*i)(* - < c(t - tJ-iQofa)- 1 
= cQ 1 (t 1 y 1 £ r (t-t 1 )- 1 . 

since Q (ti) = £oQi{ti). 
So, 

[giW 1+r giW] fc <o(^)(i-*i)- fe - 

Hence, for k > 1, or for m > 2 + r + 1 = r + 3, 



[ QoQ?- l <0{E ) 
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so 



lfo(«o)ft 2m | <O(£ )Q Ql 
If * > t , t < t u Q (t) < ££Qi(t); therefore 

f t QoQr'ds = J t T Q QiQr 3 ds = (^Qo(s) + Ro(s)) Qr 3 (s) 

-f(^ + />^) gr4(m - 3) ^- 



The first term on the right hand side is local and high order. The second term on 
the right hand side is bounded by 

c t QlQ?~ 2 ds + [ T (•••) + h.o.t. < O(S )Qj + 0(SS) {ty 3 + c t Q 2 Q?- 2 ds 

Jt Jt x Jt 

< 0(Sq) [Ql + {ty 3 ] + cS r J' 1 Q QT' l ds. 

Hence, 



T 



QoQids < O(S ) Qi + (t}~ 3 + O(8o) / Q Q? ds 



which implies that 



[ QvQ n r l ds<o{£ Q ) [Qi + (ty 



(12.63) 



(12.64) 



for all t > t . 

For < t < t , we need to estimate 



2 al(t) J' O(a 2 )p 2m ds. 



Using that for t < t ,\a (t)\ 2 < k (t) 1 < £o(t ) 1 (t) 1 the above expression is 
bounded by 

O(£ )k (t )(t)-\\n^)k (t ) < OiSM^it)- 1 (|| ln |y) (t)' 1 

<O(£ )k (t )(t)- 2 . 
This completes the proof of Proposition 12.1. 
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13. Ri\r)b) terms of Proposition 11.1 

Proposition 13.1. Assume either t < to or Montonicity Property Q on [to, to + 
5*]. Then, the terms Ri\n b ], i = 0,1 in (11.1) satisfy the estimates: 

R?[Vb] | < O(£o) [2TQoQ\ + (* > _1 (t>- 2 (Poly [P (0), Pi(0), P (T), P^T)}) 

(13.1) 

where Poly [• ■ •] stands for polynomial in the bracketed variables. 

Proof: The contributions to Ri\r}b) comes form linear and nonlinear terms 
in rjb in the Pi equations. 

Consider first the nonlinear contributions: 
In the P equations we have terms like (7.41) 

e iA +WMx,M 2 >, (xAVb\\)a . (13.2) 

Since for t G Iq, l-i we have time decay of either a or (5\ respectively, the main 
contribution is when t e I\ — [to^i]- 
For t e h the bound on r/b we have is 



\Vb 



H k <cQl /2 (t)j^-^QT(s). 



(t 

The second and third terms can be integrated by parts leaving terms of the type 

0{(5\) 0(ao)dT{xM 2 + vD- 

We also need to integrate by parts the two other terms. For this we need to pull 
out a phase factor from the leading nonlocal. 

Pulling a phase as in the proof of Proposition 12.1 we are left with estimating 
term of the type 

0{al + aoPJix^bdfrb) + (x, \Vb\ 2 dfn b )O(a ). 

As in the estimates of Proposition 10.3 , for t e h, 

ds 



\\d^ b \\ HS <CQl /2 (t)J o j^QT(s). 



with k' large for k large. 

To this end we use the following. 
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Proposition 13.2. Fort e h: 

« ds 



uQT(s) < O(£ )Ql(t). 



ho (t- S } k ' 

Proof of Proposition 13.2: 

The second term on the right hand side can be estimated, using that Q\ = 

i 

— rQoQi + CKQoQ™) an( i monotonicity of Qo, as follows: (s < £ < t) 



Repeating this argument, we have 



< o(f„)Q5(t). 

for k' — 2j > 1, which proves the Proposition 13.2. 
This Proposition together with the estimate 

H\\w^ < O(S )Ql /2 (13.3) 

implies that 

< O(S ) [ QoQl + higher order terms ] . (13.4) 

The estimates of are similar. 

It remains to estimate the linear rji, terms in the Pq, Pi eqs. 

The leading order source term of was estimated in Proposition 12.1. It 
remains to estimate the higher order corrections. 



i 
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To this end we need to estimate terms of the following type appearing in the 
Pq equation, (11.1): 

«o (\Pi\ 2 + aoPi) (xj'e-^-'^Pct^Vbds ) and 

(|A| 2 a + a a f3 1 ) (x, J* e'™ ^ P c \r) b \\ds } , 

and similar terms in the Pi equation. 

Again we focus on t E I±. Since for < s < t a and rj b are of order , 
clearly these contributions are of order 

O(£ ) [Q Ql + (t )- 3 

so it remains to estimate the s-integrals above on Ii. 

But on 7i, < 0(£o)Qo(t) < O(£ )Qi(t) and since Q is monotonic in- 
creasing on Ii, the above nonlocal term are bounded by 

0(£o)Q$(t). 

So, 

O(\Pi\ 2 a )O(£ )Ql(t) < O(£ )QlQ 
0(\a \ 2 Pi)0(£o)Ql(t) < O(£ )QlQ 
since Q < O(£ )Qi on I lm 

14. Bootstrapping it all 

We assume that to < oo, where t is given by (8.5). Consider the equations for P 
and Pi, (11.1) displayed in Proposition 11.1. Explicit in (11.1) are terms which 

(1) are driven by the dispersive part of the initial data: -Robto] and -Rq[^o] 

(2) encompass interactions of the two bound states and dispersive waves: Ri[r)b] 
and 

(3) encompass interactions between bound states: R^Pq, Pi]. 
By Proposition 11.2 the R J [r) ], J = 0, 1 terms satisfy: 

Ri[ m ] = OfattMx^it)- 2 + O(£ )2TP P? (14.1) 
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Therefore, by Proposition 9.2 and its proof (see eqn. (9.22) ), it is natural to 
introduce the functions: 

Qo(t) = Po-|y, Qi(t) = Pi + |y (14.2) 

where 

k = b + O(S )b\ , h = 10&i (14.3) 

b = (to)- 1 ( Mx + cXl ) , (14.4) 

6i = (to)- 1 ' ( MI + dXl ) , (14.5) 

where the constants c* and d* are to be chosen. We find, for any m > 4 and all 
t > 0: 

dQ ° > 2TQ Ql + R { #[r ]b ] + R^*[Q ,Q 1 ] + ^ 



dt ~ TOi ' 1 L,0J ' 2 Ln;u '^ iJ ' (t ) (t) 2 



r < -4rg Q? + ^ # M + i?2' # [Qo,Qi] + ^o;m)Jgo<5r - 



(14.6) 

The analogous reverse inequalities hold as well with slightly different constants. 

By the definition of to, Qo{to) > 0. Furthermore, using the energy estimate on 
the bound state amplitudes and (14.2) of section 5.2, we have 

Qo{t) + Qi(t) <CS , t>0. (14.7) 



We now introduce a set of norms. The norm of q(t) = ( Qo(t),Qi(t),rjb(t) ) is 
defined as 

\\q(t)\W = \Qo\vow + \Qi\ yi (t) + \\vb\\y 2 (t)- (14.8) 

The norm, ||<7(£)||y) encodes all the estimates for Qo,Qi and i] b in the intervals 
I , ii and I 2 through the following: 

|QoL(t) = SU P lQo(s)| + sup (* )(s)|Qo(s)| (14.9) 

0<s<t 0<s<min{t,t } 

\Qi\ yi{t) = sup |Qi(s)| + sup |s-ti|Qi(s)r'Q (ti)Qi(ti) (14.10) 

0<s<t ti<s<t 

\\Vb\\y 2 (t) = SU P (s)ll^(s)llw*.<» + SU P \\Vb(s)\\ W k,^ 

0<s<min{i,io} to<s<min{t,ti} 

+ sup (s — ti)* \\r] b (s)\\ W k,oo (14-11) 

t!<S<t 
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In these definitions we use the convention that terms for which the s- range empty 
are set to zero. 



By the H 1 a priori bounds 

sup |Q (s)| < So (l +£olM|oc) (14.12) 



0<s<t 



and by definition of t , (8.5), for t < oo, 

sup (t )(s)\Q (s)\ <Us + Mx ) • (14.13) 

0<s<to * 

In terms of these norms, we have bounds on Rf*. By Proposition 13.1 
h -[Vb\\ ^ ^6;||gwi|yVoWi-i-o 

By Proposition 12.1 



R Av b }\ + \Rl'*h]\ < o(sS)\\qm^QoQl + c ^l\)!^ . (14.14) 

WW 



\R° 2 '*[Qo,Qi]\ < O(E§)Q Ql + g ^tfff* (14-15) 
l^2' # [Qo,Qi]| < O(££)Q Ql + C ^ + }^}f" (14.16) 



where / > 2. 



By the definition of Iq, (0 < t < t ) and Propositions 11.2, 12.1 and 13.1, we 
have estimates (14.14,14.15,14.16). Therefore, for an appropriate choice of c* and 
we have for < t < t 

— > 2 rg g 1 + 

< -4r'Q Q? + O(^o;m) v /^gr - #-^7tt, (14.17) 



where m > 4. Note that by definition of J , Qo(^) < for t e 7 - 

By continuity, (14.17) holds for t < t < t + 5, for some 5 > 0. It follows, using 
that Qo(t ) > and Propositions and , that (11.2) (Monotonicity Property 
Q) holds on t < t < t + 5. Therefore, by Propositions 11.2, 12.1 and 13.1 the 
terms J and Ji in (9.8) both satisfy the bound 

\Jk\ <ttT77a/ + ZoQoQl, p>0,k = 0,l- (14.18) 

WW 
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Therefore, for S sufficiently small (14.17) holds with c*, replaced by c*/2, d*/2. 
Define 

T* = sup{t>t : (14.17) holds for some c* > and d* > 0} (14.19) 

For t E [0, T*), ||?(0I|y is small. We claim that T* = oo. Suppose t < ^* < °°- 
Then, for £ G [£o,T*) we have, by Propositions 14 and 14, that the monotonicity 
property (11.2) holds at t = T* and slightly beyond. Thus, the a priori bounds 
on Jo, Ji of Propositions 11.2, 12.1 and 13.1, the rjb - bounds of Propositions 10.7, 
10.8 and 10.9, and the smallness of Qo and Q\ imply persistence of the inequalities 
(14.17), with perhaps a slightly smaller choice of positive constants c* and d*. This 
implies that T* = oo. 

15. Nongeneric behavior 

Recall that t is defined by (8.5). and consider the case where t — oo . We would 
like to show that 

Po(t) as t — > oo 

Pi(t) has a limit. 

The following is a consequence of the definition of to- 

Proposition 15.1. Assume to = oo. Then, Po(t) = O ([i]o]x(t)~ 2 ) ■ Therefore, 
ao — * and the ground state decays. 



Proposition 15.2. Let t = oo. Then, (3\ has a limit as t — > oo. 

States with this (nongeneric) behavior were constructed in [50] . 
Proof of Proposition 15.2: The equation (7.7), together with the above esti- 
mate 

p (t) = o(m)(t)- 2 (i5.i) 

implies 

dP 1 



dt 



= {-ATPoPl + Om(t)- 3 / 2 )(l + O(P ) + O(Pi)) 
+ 3?(ce iA+t |A| 2 Aao) + h.o.t. (15.2) 



97 



To show that Pi has a limit, we show that /q d s P2(s)ds has a limit. All terms 
other than the O(ao) term, on the right hand side are absolutely integrable since 
Po = 0((t)-*). 

It is left to integrate the Oflftp/^ao) term. For T given, let fi\ = (3i(T) 2 . 
Then, equation (7.52) reads 

2id t ao = A(0o*, ^i{t) 2 )ao + integrable in t. (15.3) 

Using the expression for \J>! (t) = aitpi(-, |«i| 2 ): 

2id t a = A(^o*, ai(t) 2 tpl # )a + h.o.t. 

= A(^o*, Oe- 2lA+ ^i(t) 2 « (t) + 0(a 2 (T)^) + h.o.t. 

= Xe-^^it) + Ae^+^W - P 2 (T)]a (t) + 0(±) 

+ h.o.t. (15.4) 

where A = A^V'L) and = (3f(T). 

Solving the homogeneous part of (15.4): 

2id t a = Xe-^+'pl&it) (15.5) 

we have, using the Ansatz 

a (t) = A(t)e l{x ~ a)t + B(t)e l{x+a)t 

with 

A ~ Ae- 2iA+ *[/? 2 (t) - (3 2 (T)]Ae i6 ^ + h.o.t.. 

and a similar equation for P(£). 
We have 
dP 

— 1 = -4rP P 2 + 3? (ce^'^il^A + e^^il^P ) + h.o.t., 

6 U, 6b 7^ 0. Integration of the above equation, integration by parts (twice) of the 
A and B, implies 

(t) 1 *\p 1 (t)-p 1 (T)\ <c{t)t jf T | ) 9 1 | !, ]9 1 (0- 1 <0[^(0-^(r)raO* , + h. .t.<t>* 

< C£ T n ( [ T (ty^dA ( sup (i3)|P 1 (i / )-P 1 (r)A +SU P (t)5 h.o.t. 

\Jt J \0<t'<T J t 

^\P 1 (t)-P 1 (T)\<C(t)- 1 * 
which implies integrability of A(t) and limit of Pi(t). 
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16. Notation 



3tz and ^sz denote the real and imaginary parts of a complex number z 
z denotes the complex conjugate of z 



( x ) = y/l + \x\ 2 , t > 

P Cif - projection onto the continous spectral part of the self-adjoint operator, H. 

(f,g) =S7 9 



a \ _ l a 
c.c. I [a 



;i6.i) 



For j = 1, 2, let ttj : C 2 -> C 1 be defined by: 
Plemelj identities 

(x^iO)" 1 = P.V. x' 1 ± i7r5(a;) (16.3) 

x(x,p) denotes are real- valued localized function of x which depends smoothly on 
a parameter, p. 

xj^ denotes a spatially localized function of order \aj\ k , as — > 0. 

Oj^ denotes a quantity which is of order \aj\ k as jctj | -> 0. Both Xk ] and Oj? are 

invariant under the map otj i— > aje n . 

or^Offoff, fc = * 1 + k. 

17. Appendix B - Proof of Proposition 4.2 

Parts (1), (3) and (5) of Proposition 4.2 follow from [54]. We now prove parts (2) 
and (4) by a perturbation argument about the case cto = 0. 

— * — * 

Consider the eigenvalue problem Tl f = {if. Since a is assumed small it is 
natural to make explicit the leading order and perturbation terms. Thus we have 



(H-E *)-EU\a \*I + ti( 2 £f ^ 



Recall that -Eg 1 '* and are defined in Proposition 3.1. 



/ = (17.1) 
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The zeroth order problem (a = 0) is 

a 3 (#-£ *)/o = /Wo, (17.2) 
which has two linearly independent solutions: 

Ho = S*i - S *, fo = ( J J ipu (17.3) 



/i = -(£7.i -So*), o-i /o- (17-4) 

We develop the perturbation theory of (17.3). That of the second is completely 
analogous. 

For a and small we define the perturbations about the zeroth order eigenstates 

via: 

/ = fo + fi (17.5) 

li = Eu-Eot + fa (17.6) 

Substitution into (17.1) yields: 

[a 3 (H - Eq.) - (E u - Eq,)I ]fr 

+ |ao| 2 S$Wi - \^ a 3 ( °* ) /i + //i/i (17.7) 

We consider, individually, the first and second equations of the system (17.7), 

— * 

governing /y = ftjfi, j — 1, 2. The first component of (17.7) is: 

(H-E u )fu = |«o| 2 (si 1} -2A^ 2 )^i* + /ii^i* 

+ |ao| 2 (^ 1) -2A^ 2 )/ 11 -Aa^2/ 12 + /(il / 11 (17.8) 

Let z/* = 2S * — Si*. The second component of (17.7) is: 



(#-0/i2 = -Ac^VoVi* + Kl 2 ^/: 



12 



-A^ 2 (c^ 2 /ii + 2|ao| 2 /i 2 ) - fiifis (17.9) 
We wish to make the dependence of /i on a and TZq explicit. Define 

Hi = \a \ 2 fli, fn = |ao| 2 /n, /12 = Oo 2 /i2 (17.10) 
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Equations (17.8) and (17.9) reduce to the following system for /n and fi 2 : 

(H-E u )fu = (e^ - 2A^ 2 ) Vi* + AiV'i* 

+ |« | 2 - 2A^ 2 ) / n - A|a | 2 ^o/i2 + KlWii(17.11) 
(H-v*)f 12 = -^u + \a \ 2 E^f 12 

-A|« |Vo (2/ii + 2/12) - Kl W12 (17.12) 

We seek a solution to the system (17.11), (17.12): 

Kl 2 ^ (/ii(|a o r),/i2(|ao| 2 ),/i(l«o| 2 )) G L 2 x L 2 x IR (17.13) 

defined in a neighborhood of «o — 0. 

For a = the system (17.11), (17.12) reduces to: 

(H-Eu)^ = (4 1) (0)-2A^)^i*+/ii^ (17.14) 
(H-v,)f° 2 = -A^i* (17.15) 

Note that z/* < is not in the spectrum of H. Therefore, (17.15) is solvable for 
fi2 and we have: 

/° 2 = -X(H - ^rVo*^*- (17.16) 

Since (if — E u )ipu = 0, (17.14) is solvable if and only if its right hand side is 
orthogonal to ipi*. This determines /ii(0): 

h(0) = -^ 1} (0) + 2A(^,^>, (17.17) 

and now (17.14) can be solved for f^. 

To solve in a neighborhood of a = we proceed as follows. Rewrite (17.12) 

as 

(H + K| 2 ^i2 - v*)fu = -A^oVi* - A|a | Wn, (17.18) 
where W i2 is a multiplication operator defined by 

W 12 (\a \\jli) = -E^ +2A^ 2 +h- (17-19) 

For jli in a fixed compact set and «o sufficiently small, the operator H+ \ao\ 2 Wi 2 — 
has a bounded inverse, -B(|o;o| 2 ). Thus, 

/l2 = /l2[/ll, |«oH 

= -AS(|a | 2 )V'o^i.-A|ao| 2 S(|aor)^o/ii (17-20) 
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Substitution of (17.20) into (17.11) yields the following closed equation for fa: 

(H - E u )fa = (£<j 1} - 2A^ 2 ) Vi. + Mi* + \u \ 2 W n fa, (17.21) 

where the operator Wu is defined by 

W u (\a \ 2 ,fa) = (E™ - 2A^ 2 ) + fa + Vo/i2[-, |«o| 2 ]- (17.22) 

Setting the inner product of the right hand side of (17.21) equal to zero, gives the 
solvability condition for (17.21): 

fa = 2\(i> 2 ,i>l) - - \a \ 2 (ij u , W n fa) (17.23) 

The system (17.21), (17.23) is of the form: 

F(fa,fa,s) = 0, (17.24) 

with the solution fa = fa, fa = fa(0),s = defined by (17.17). Furthermore, 
the Jacobian of F(fa,fa, s) with respect to (fa, fa) evaluated at (fa,fa(Q), 0) is 
given by: 



V 1 ) 

which maps H 2 x 1R one to one and onto L 2 x {(g, : g G L 2 }. Therefore, 
by the implicit function theorem [31], we have a real analytic curve of solutions 
s i — y (fa(s),fa(s),s), defined in a neighborhood of s = and coinciding with 
(fii,fa(0), 0) for s = 0. The family of solutions we seek is obtained by restriction 
to s = \a \ 2 > 0. This completes the proof of Proposition 4.2. 

18. Appendix C: A commutator term 

In this section we record a calculation of a "commutator term" appearing in the 
modulation equations of section 5. 

Proposition 18.1. 

i<$(<7 3 £oi),$2) = id t (\a \ 2 ) (f o Go[ \ ) > ®i) 

+ d ao \a \ 2 (xG ^ 1 ) ' $2 ) (18 - x) 

i<$(<73&2),$2) = ^(Kl 2 ) ^W?o( \ ) , 

+ (dao) |« | 2 (x^Go ( \ J , $ 2 ) (18.2) 
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proof: By direct computation from (4.23) 

d t G (t)=i(d tlo )(T3Go(t). (18.3) 

Note also that by (4.37) 

0-3G0 ( 1 ) F ° = ^ 01 = 2< ^ 01 + \ a o\ 2(T 3 G o ^ J j X(x; \a \ 2 ) 
G ^ j j = £02 = ^Co2 + K| 2 G ^ J x(a:; K| 2 ) 
Using these relations we have for j = 1 that 

$(<7 3 £oi(*)) = ^(l«o| 2 )G 
+ i(9 t 7o)cr 3 Go 
= ft(K| 2 )G 

Substitution into the inner product (<9 t (o" 3 £oi(£)), $2) and using the constraint 
(Coi(£), $2) = yields the result for j = 1. 
For j = 2 

$(<7 3 $)2(f)) = ^(K|> 3 G ( J ^o(|«o| 2 ) 

+ ^(%o)G ( j ) i^(K| 2 ) 

= d t (\a \ 2 )a 3 G ( j j ^o(K| 2 ) + ^(%o)Co 2 (t) + %o|«o| 2 x(^; KIV3G 

Substitution into the inner product (<9t(o" 3 £o2(£)), $2) and using the constraint 
(Co2(t), $2) = yields the result for j = 2. This completes the proof of Proposition 
18.1. 
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